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The LIGO-II gravitational-wave interferometers (ca. 2006-2008) are designed to have sensitivities 
near the standard quantum limit (SQL) in the vicinity of 100 Hz. This paper describes and analyzes 
possible designs for subsequent, LIGO-III interferometers that can beat the SQL. These designs 
are identical to a conventional broad-band interferometer (without signal recycling), except for 
new input and/or output optics. Three designs are analyzed: (i) a squeezed-input interferometer 
(conceived by Unruh based on earlier work of Caves) in which squeezed vacuum with frequency- 
dependent (FD) squeeze angle is injected into the interferometer's dark port; (ii) a variational- 
output interferometer (conceived in a different form by Vyatchanin, Matsko and Zubova), in which 
homodyne detection with FD homodyne phase is performed on the output light; and (iii) a squeezed- 
variational interferometer with squeezed input and FD-homodyne output. It is shown that the FD 
squeezed-input light can be produced by sending ordinary squeezed light through two successive 
Fabry-Perot filter cavities before injection into the interferometer, and FD-homodyne detection 
can be achieved by sending the output light through two filter cavities before ordinary homodyne 
detection. With anticipated technology (power squeeze factor e~ 2R — 0.1 for input squeezed vacuum 
and net fractional loss of signal power in arm cavities and output optical train e„ = 0.01) and using 
an input laser power I in units of that required to reach the SQL (the planned LIGO-II power, 
^sql), the three types of interferometer could beat the amplitude SQL at 100 Hz by the following 
amounts fi = y/Sh/ \J and with the following corresponding increase V = l//i 3 in the volume 
of the universe that can be searched for a given non-cosmological source: 
Squeezed-Input — fi ~ V e~ 2R ~ 0.3 and V ~ 1/0. 3 3 ~ 30 using / /7sql = L 

Variational- Output — fj, ~ e^/ 4 ~ 0.3 and V ~ 30 but only if the optics can handle a ten times 
larger power: 7 //sql — l/s/u = 10. 

Squeezed- Varational — \i = 1.3(e -2ii e„) 1 ' /4 ~ 0.24 and V ~ 80 using Lj/^sql = 1; and /i ~ 
(e- 2S e,) 1/4 ~ 0.18 and V ~ 180 using J // S ql = y/e- 2R /e* ~ 3.2. 

04.80. Nn, 95.55.Ym, 42.50.Dv, 03.65.Bz 
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I. INTRODUCTION AND SUMMARY 

In an interferometric gravitational- wave detector, laser 
light is used to monitor the motions of mirror-endowed 
test masses, which are driven by gravitational waves h(t). 
The light produces two types of noise: photon shot noise, 
which it superposes on the interferometer's output sig- 
nal, and fluctuating radiation-pressure noise, by which 
it pushes the test masses in random a manner that can 
mask their gravity-wave-induced motion. The shot-noise 
spectral density scales with the light power I entering 
the interferometer as Sf^ ot oc 1/I ; the radiation-pressure 
noise scales as S 1 ^ oc I . 

In the first generation of kilometer-scale interferome- 
ters (e.g., LIGO-I, 2002-2003 §)), the laser power will 
be low enough that shot-noise dominates and radiation- 
pressure noise is unimportant. Tentative plans for the 
next generation interferometers (LIGO-II, ca. 2006-2008) 
include increasing I a to the point that, = Sf^ ot 
at the interferometers' optimal gravitational-wave fre- 
quency, 0/27T ~ 100 Hz. The resulting net noise Sh = 

S h P + ss h 0t = 2SS h 0t is the lowest that can be achieved 
with conventional interferometer designs. Further in- 
creases of light power will drive the radiation-pressure 
on upward, increasing the net noise, while reductions of 
light power will drive the shot noise upward, also increas- 
ing the net noise. 

This minimum achievable noise is called the "Standard 
Quantum Limit" (SQL) g and is denoted fif QL = /i§ QL . 
It can be regarded as arising from the effort of the quan- 
tum properties of the light to enforce the Heisenberg 
uncertainty principle on the interferometer test masses, 
in just the manner of the Heisenberg microscope. In- 
deed, a common derivation of the SQL is based on the 
uncertainty principle for the test masses' position and 
momentum ||]: The light makes a sequence of measure- 
ments of the difference x of test-mass positions. If a 
measurement is too accurate, then by state reduction it 
will narrow the test-mass wave function so tightly (Ax 
very small) that the momentum becomes highly uncer- 
tain (large Ap), producing a wave- function spreading 
that is so rapid as to create great position uncertainty 
at the time of the next measurement. There is an op- 
timal accuracy for the first measurement — an accuracy 
that produces only a factor \[2 spreading and results in 
optimal predictability for the next measurement. This 
optimal accuracy corresponds to hsQh- 

Despite this apparent intimate connection of the SQL 
to test-mass quantization, it turns out that the test-mass 
quantization has no influence whatsoever on the output 
noise in gravitational- wave interferometers Q] . The sole 
forms of quantum noise in the output are photon shot 
noise and photon radiation-pressure noise.F] 



1 In brief, the reasons for this are the following: The inter- 



Vladimir Braginsky (the person who first recognized 
the existence of the SQL for gravitational-wave detectors 
and other high-precision measuring devices Q) realized, 
in the mid 1970s, that the SQL can be overcome, but to 
do so would require significant modifications of the ex- 
perimental design. Braginsky gave the name Quantum 
Nondemolition (QND) to devices that can beat the SQL; 
this name indicates the ability of QND devices to pre- 
vent their own quantum properties from demolishing the 
information one is trying to extract Q . 

The LIGO-I interferometers are now being assembled 
at the LIGO sites, in preparation for the first LIGO 
gravitational-wave searches. In parallel, the LIGO Sci- 
entific Community (LSC) is deeply immersed in R&D 
for the LIGO-II interferometers [[7j, and a small portion 
of the LSC is attempting to invent practical designs for 
the third generation of interferometers, LIGO-III. This 
paper is a contribution to the LIGO-III design effort. 

In going from LIGO-II to LIGO-III, a large number of 
noise sources must be reduced. Perhaps the most serious 
are the photon shot noise and radiation pressure noise 
("optical noise"), and thermal noise in the test masses 
and their suspensions ^ n this paper we shall deal 

solely with the shot noise and radiation pressure noise 
(and the associated SQL); we shall tacitly assume that 
all other noise sources, including thermal noise, can be 
reduced sufficiently to take full advantage of the optical 
techniques that we propose and analyze. 

Because LIGO-II is designed to operate at the SQL, 
in moving to LIGO-III there are just two ways to reduce 
the optical noise: increase the masses m of the mirrored 
test masses (it turns out that hgq L oc 1/m), or redesign 
the interferometers so they can perform QND. The tran- 
sition from LIGO-I to LIGO-II will already (probably) 



ferometer's measured output, in general, is one quadrature of 
the electric field [the b( of Eqs. ( p4| ) and jlo| ) below] , and this 
output observable commutes with itself at different times by 
virtue of Eqs. (Q) with a — > b. This means that the digi- 
tized data points (collected at a rate of 20 kHz) are mutually 
commuting Hermitian observables. One consequence of this 
is that reduction of the state of the interferometer due to 
data collected at one moment of time will not influence the 
data collected at any later moment of time. Another conse- 
quence is that, when one Fourier analyzes the interferometer 
output, one puts all information about the initial states of the 
test masses into data points near zero frequency, and when 
one then filters the output to remove low-frequency noise 
(noise at / = Q./2n < 10 Hz), one thereby removes from 
the data all information about the test-mass initial states; 
the only remaining test-mass information is that associated 
with Heisenberg-picture changes of the test-mass positions at 
/ > 10 Hz, changes induced by external forces: light pressure 
(which is quantized) and thermal- and seismic-noise forces 
(for which quantum effects are unimportant). See Ref. for 
further detail. 
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entail a mass increase, from m — 11 kg to m = 30 kg, in 
large measure because the SQL at 11 kg was unhappily 
constraining [Q. Any large further mass increase would 
entail great danger of unacceptably large noise due to en- 
ergy coupling through the test-mass suspensions and into 
or from the overhead supports (the seismic isolation sys- 
tem); a larger mass would also entail practical problems 
due to the increased test-mass dimensions. Accordingly, 
there is strong motivation for trying to pursue the QND 
route. 

Our Caltech and Moscow-University research groups 
are jointly exploring three approaches to QND interfer- 
ometer design: 

• The conversion of conventional interferometers into 
QND interferometers by modifying their input 
and/or output optics [this paper]. This approach 
achieves QND by creating and manipulating cor- 
relations between photon shot noise and radiation 
pressure noise; see below. It is the simplest of our 
three approaches, but has one serious drawback: 
an uncomfortably high light power, W c - lrc > 1 MW, 
that must circulate inside the interferometers' arm 
cavities M. It is not clear whether the test- mass 
mirrors can be improved sufficiently to handle this 
high a power in a sufficiently noise-free way. 

• A modification of the interferometer design (includ- 
ing using two optical cavities in each arm) so as to 
make its output signal be proportional to the rela- 
tive speeds of the test masses rather than their rel- 
ative positions fl(],|ll| . Since the test-mass speed is 
proportional to momentum, and momentum (un- 
like position) is very nearly conserved under free 
test-mass evolution on gravity-wave timescales (~ 
0.01 sec), the relative speed is very nearly a "QND 
observable" 0] and thus is beautifully suited to 
QND measurements. Unfortunately, the resulting 
speed-meter interferometer, like our input-output- 
modified interferometers, suffers from a high circu- 
lating light power ||, W circ > 1 MW. 

• Radical redesigns of the interferometer aimed at 
achieving QND performance with W C1IC well below 
1 MW These, as currently conceived by Bra- 
ginsky, Gorodetsky and Khalili, entail transfering 
the gravitational-wave signal to a single, small test 
mass via light pressure, and using a local QND sen- 
sor to read out the test mass's motions relative to 
a local inertial frame. 

In this paper we explore the first approach. The foun- 
dation for this approach is the realization that: (i) photon 
shot noise and radiation-pressure noise together enforce 
the SQL only if they are uncorrelated; see, e.g., Ref. [Q; 
(ii) whenever carrier light with side bands reflects off a 
mirror (in our case, the mirrors of an interferometer's 
arm cavities), the reflection ponderomotively squeezes the 
light's side bands, thereby creating correlations between 
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FIG. 1. Schematic diagram of a squeezed-input interfer- 
ometer. 



their radiation-pressure noise in one quadrature and shot 
noise in the other; (iii) these correlations are not ac- 
cessed by a conventional interferometer because of the 
particular quadrature that its photodiode measures; (iv) 
however, these correlations can be accessed by (concep- 
tually) simple modifications of the interferometer's input 
and/or output optics, and by doing so one can beat the 
SQL. These correlations were first noticed explicitly by 
Unruh fl4| , but were present implicitly in Braginsky's 
earlier identification of the phenomenon of ponderomo- 
tive squeezing JL5|Jl6| . 

In this paper we study three variants of QND interfer- 
ometers that rely on ponderomotive-squeeze correlations: 

(i) Squeezed- Input Interferometer: Unruh fl4|| (build- 
ing on earlier work of Caves ]l7|]) invented this design 
nearly 20 years ago, and since then it has been reana- 
lyzed by several other researchers Jl^,|l9|. In this design, 
squeezed vacuum is sent into the dark port of the in- 
terferometer ( "modified input" ) and the output light is 
monitored with a photodetector as in conventional inter- 
ferometers. 

For a broad-band squeezed-input interferometer, the 
squeeze angle must be a specified function of frequency 
that changes significantly across the interferometer's op- 
erating gravity-wave band. (This contrasts with past 
experiments employing squeezed light to enhance inter- 
ferometry |2^j2l[], where the squeeze angle was constant 
across the operating band.) Previous papers on squeezed- 
input interferometers have ignored the issue of how, 
in practice, one might achieve the required frequency- 
dependent (FD) squeeze angle. In Sec. VC, we show that 
it can be produced via ordinary, frequency-independent 
squeezing (e.g., by nonlinear optics p2|), followed by fil- 
tration through two Fabry-Perot cavities with suitably 
adjusted bandwidths and resonant-frequency offsets from 
the light's carrier frequency. A schematic diagram of the 
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FIG. 2. Schematic diagram of a squeezed- variational in- 
terferometer. A variational-output interferometer differs from 
this solely by replacing the input squeezed vacuum by ordi- 
nary vacuum. 



resulting squeezed-input interferometer is shown in Fig. 
[I] and is discussed in detail below. Our predicted per- 
formance for such an interferometer agrees with that of 
previous research. 

(ii) Variational- Output Interferometer: Vyatchanin, 
Matsko and Zubova invented this design conceptually in 
the early 1990's [p3| p5f. It entails a conventional inter- 
ferometer input (ordinary vacuum into the dark port), 
but a modified output: instead of photodetection, one 
performs homodyne detection with a homodyne phase 
that depends on frequency in essentially the same way 
as the squeeze angle of a squeezed-input interferometer. 
Vyatchanin, Matsko and Zubova did not know how to 
achieve FD homodyne detection in practice, so they pro- 
posed approximating it by homodyne detection with a 
time-dependent (TD) homodyne phase. Such TD ho- 
modyne detection can beat the SQL, but (by contrast 
with FD homodyne) it is not well-suited to gravitational- 
wave searches, where little is known in advance about the 
gravitational waveforms or their arrival times. In this 
paper (Sec. ^ and Appendix |^), we show that the de- 
sired FD homodyne detection can be achieved by sending 
the interferometer's output light through two successive 
Fabry-Perot cavities that are essentially identical to those 
needed in our variant of a squeezed-input interferometer, 
and by then performing conventional homodyne detec- 
tion with fixed homodyne angle. A schematic diagram of 
the resulting variational-output interferometer is shown 
in Fig. |. 

(Hi) Squeezed- Variational Interferometer: This design 



(not considered in the previous literature^]) is the obvious 
combination of the first two; one puts squeezed vacuum 
into the dark port and performs FD homodyne detection 
on the output light. The optimal performance is achieved 
by squeezing the input at a fixed (frequency-independent) 
angle; filtration cavities are needed only at the output 
(for the FD homodyne detection) and not at the input; 
cf. Fig. | 



In Sec. IV we compute the spectral density of the 
noise for all three designs, ignoring the effects of opti- 
cal losses. We find (in agreement with previous anal- 
yses [[18 19 1) that, when the FD squeeze angle is op- 
timized, the squeezed-input interferometer has its shot 
noise and radiation-pressure noise both reduced in am- 
plitude (at fixed light power) by e~ R , where R is the 
(frequency-independent) squeeze factor; see Fig. 2 be- 
low. This enables a lossless squeezed-input interferom- 
eter to beat the SQL by a factor e~ R (when the power 
is optimized) but no more. By contrast, the lossless, 
variational-output interferometer, with optimized FD ho- 
modyne phase, can have its radiation-pressure noise com- 
pletely removed from the output signal, and its shot noise 
will scale with light power as \j\[T as for a conven- 
tional interferometer. As a result, the lossless variational- 
output interferometer can beat the SQL in amplitude by 
\/^sql/2/o, where /sql is the light power required by a 
conventional interferometer to reach the SQL. The op- 
timized, lossless, squeezed- variational interferometer has 
its radiation-pressure noise completely removed, and its 
shot noise reduced by e~ R , so it can beat the SQL in 
amplitude by e~ R y / I S Q L /2I . 

Imperfections in squeezing, in the filter cavities, and 
in the homodyne local-oscillator phase will produce er- 
rors AA in the FD squeeze angle A(fi) of a squeezed- 
input or squeezed- variational interferometer, and A£ in 
the FD homodyne phase ((£Y) of a variational-output or 
sq u eezed- variational interferometer. At the end of Sec. 



VIE , we shall show that, to keep these errors from seri- 
ously compromising the most promising interferometer's 
performance, |AA| must be no larger than ~ 0.05 radian, 
and |A£| must be no larger than ~ 0.01 radian. This 
translates into constraints of order five percent on the 
accuracies of the filter cavity finesses and about 0.01 on 
their fractional frequency offsets and on the homodyne 
detector's local-oscillator phase. 

The performance will be seriously constrained by un- 
squeezed vacuum that leaks into the interferometer's op- 
tical train at all locations where there are optical losses, 



2 A design similar to it has previously been proposed and 
analyzed |w| for a simple optical meter, in which the position 
of a movable mirror (test mass) is monitored by measuring 
the phase or some other quadrature of a light wave reflected 
from the mirror. In this case it was shown that the SQL can 
be beat by a combination of phase-squeezed input light and 
TD homodyne detection. 
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whether those losses are fundamentally irreversible (e.g. 
absorption) or reversible (e.g. finite transmissivity of an 
arm cavity's end mirror). We explore the effects of such 



optical losses in Sec. VI. The dominant losses and as- 
sociated noise production occur in the interferometer's 
arm cavities and FD filter cavities. The filter cavities' 
net losses and noise will dominate unless the number of 
bounces the light makes in them is minimized by mak- 
ing them roughly as long as the arm cavities. This sug- 
gests that they be 4km long and reside in the beam tubes 
alongside the interferometer's arm cavities. To separate 
the filters' inputs and outputs, they might best be trian- 
gular cavities with two mirrors at the corner station and 
one in the end station. 

Our loss calculations reveal the following: 
The squeezed-input interferometer is little affected by 
losses in the interferometer's arm cavities or in the out- 
put optical train. However, losses in the input optical 
train (most seriously the filter cavities and a circulator) 
influence the noise by constraining the net squeeze fac- 
tor e~ 2R of the light entering the arm cavities. The re- 
sulting noise, expressed in terms of e~ 2R , is the same 
as in a lossless squeezed-input interferometer (discussed 
above): With the light power optimized so I Q = /sql, 
the squeezed-input interferometer can beat the ampli- 
tude SQL by a factor » = ^/SiJ^S^ ~ Ve- 2R ~ 0.3 

(where e~ 2R ~ 0.1 is a likely achievable value of the 
power squeeze factor). 

The variational- output and squeezed- variational inter- 
ferometers are strongly affected by losses in the inter- 
ferometer's arm cavities and in the output optical train 
(most seriously: a circulator, the two filter cavities, 
the mixing with the homodyne detector's local-oscillator 
field, and the photodiode inefficiency). The net fractional 
loss e* of signal power and (for squeezed-variational) the 
squeeze factor e~ 2R for input power together determine 
the interferometer's optimized performance: The ampli- 
tude SQL can be beat by an amount [i = (e~ 2R e t ,) 1 / 4 , 
and the input laser power required to achieve this op- 
timal performance is I /Isql — \/ e ~ 2R / e *- In par- 
ticular, the variational-output interferometer (no input 
squeezing; e~ 2R = 1), with the likely achievable loss 
level e* = 0.01, can beat the SQL by the same amount 
as our estimate for the squeezed-input interferometer, 
jj, ~ ftj^ ~ 0.3, but requires ten times higher input 
optical power, I /Isql — — 10 — which could 

be a very serious problem. By contrast, the squeezed- 
variational interferometer with the above parameters has 
an optimized performance fi ~ (0.1 x 0.01) 1 / 4 ~ 0.18 
(substantially better than squeezed-input or variational- 
output), and achieves this with an optimizing input 
power I /Isql = \/0. 1/0.01 ~ 3.2. If the input power is 
pulled down from this optimizing value to / /isQL = 1 
so it is the same as for the squeezed-input interferometer, 
then the squeezed-variational performance is debilitated 
by a factor 1.3, to fj, ~ 0.24, which is still somewhat 
better than for squeezed- input. 



It will require considerable R&D to actually achieve 
performances at the above levels, and there could be a 
number of unknown pitfalls along the way. For exam- 
ple, ponderomotive squeezing, which underlies all three 
of our QND configurations, has never yet been seen in the 
laboratory and may entail unknown technical difficulties. 

Fortunately, the technology for producing squeezed 
vacuum via nonlinear optics is rather well developed |22^] 
and has even been used to enhance the performance of 
interferometers [ poU2lf| . Moreover, much effort is being 
invested in the development of low-loss test-mass sus- 
pensions, and this gives the prospect for new (pondero- 
motive) methods of generating squeezed light that may 
perform better than traditional nonlinear optics. These 
facts, plus the fact that, in a squeezed-input configura- 
tion, the output signal is only modestly squeezed and 
thus is not nearly so delicate as the highly-squeezed 
output of an optimally performing squeezed-variational 
configuration, make us feel more confident of success 
with squeezed-input interferometers than with squeezed- 
variational ones. 

On the other hand, the technology for a squeezed- 
variational interferometer is not much different from that 
for a squeezed-input one: Both require input squeezing 
and both require filter cavities with roughly the same 
specifications; the only significant differences are the need 
for conventional, frequency-independent homodyne de- 
tection in the squeezed-variational interferometer, and 
its higher-degree of output squeezing corresponding to 
higher sensitivity. Therefore, the squeezed-variational 
interferometer may turn out to be just as practical as 
the squeezed-input, and may achieve significantly better 
overall performance at the same laser power. 

This paper is organized as follows: In Sec. || we sketch 
our mathematical description of the interferometer, in- 
cluding our use of the Caves- Schumaker ||^J2^] formalism 
for two-photon quantum optics, including light squeezing 
(cf. Appendix^); and we write down the interferometer's 
input-output relation in the absence of losses [Eq. (|l^); 
cf. Appendix [b] for derivation]. In Sec. f\% relying on 
our general lossless input-output relation (|16|), we derive 
the noise spectral density Sh(f) for a conventional inter- 
ferometer and elucidate thereby the SQL. In Sec. IV, wc 



describe mathematically our three QND interferometer 
designs and, using our lossless input-output relation ( Jlq) , 
derive their lossless noise spectral densities. In Sec. |V|, 
we show that FD homodyne detection can be achieved by 
filtration followed by conventional homodyne detection, 
and in Appendix ^| we show that the required filtration 
can be achieved by sending the light through two suc- 
cessive Fabry-Perot cavities with suitably chosen cavity 
parameters. We list and discuss the required cavity pa- 
rameters in Sec. [v| In Sec. [Vl|, we compute the effects of 
optical losses on the interferometers' noise spectral den- 
sity; our c omp utation relies on an in pu t-output relat ion 



|) and fllOlp derived in Appendix g[ In Sec. |VlI 



we 



discuss and compare the noise performan ces o f our three 
types of inteferometers. Finally, in Sec. VIII we briefly 
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FIG. 3. Gravitational-wave interferometer with two inputs 
(the carrier which has power I a entering the bright port, and 
quantum field a entering the dark port) and one relevant out- 
put (the quantum field b leaving the dark port). 



recapitulate and then list and briefly discuss a number 
of issues that need study, as foundations for possibly im- 
plementing these QND interferometers in LIGO-III. 

This paper assumes that the reader is familiar with 
modern quantum optics and its theoretical tools as pre- 
sented, for example, in Refs. p8fl. 



II. MATHEMATICAL DESCRIPTION OF THE 
INTERFEROMETER 

A. Input and Output fields 

Figure || shows the standard configuration for a 
gravitational-wave interferometer. In this subsection we 
focus on the beam splitter's input and output. In our 
equations we idealize the beam splitter as infinitesimally 
thin and write the input and output fields as functions 
of time (not time and position) at the common centers 
of the beams as they strike the splitter. 

At the beam splitter's bright port the input is a carrier 
field, presumed to be in a perfectly coherent state with 
power I a ~ 10 kW (achieved via power recycling [pj|), 
angular frequency uj ~ 1.78 x 10 15 sec _1 (1.06 micron 
light), and excitation confined to the cos(w i) quadra- 
ture (i.e., the mean field arriving at the beam splitter is 
proportional to cos(w i)). 

At the dark port the input is a (quantized) electromag- 
netic field with the positive-frequency part of the electric 
field given by the standard expression 



2nhuj 

~Ac~ 



< dui 
2^ 



(1) 



Here A is the effective cross sectional area of the beam 
and a^j is the annihilation operator, whose commutation 
relations arc 



[a u ,a ui ']=0, [a u ,a',] = 2n5(u -lj') 



(2) 



Throughout this paper we use the Heisenberg Picture, 
so E^ evolves with time as indicated. However, our 
creation and annihilation operators a u and aj, are fixed in 
time, with their usual Heisenberg-Picture time evolutions 
always factored out explicitly as in Eq. (|l|) . 

We split the field (|l|) into side bands about the carrier 
frequency w 0} u> = u> ± f2, with side-band frequencies f2 
in the gravitational-wave range ~ 60 to ~ 6000sec~ 1 (10 
to 1000 Hz), and we define 



a + 



(3) 



As in Eq. (g), we continue to use a prime on the subscript 
to denote frequency Ct': a + > = a w<J +n'- Correspondingly, 
the commutation relations (Q) imply for the only nonzero 
commutators 

[a + ,a i + ,] = 2'K5(Q-Q') , [o_, at,] = 2nd{Q - O') ; (4) 
and expression (0) for the dark -port input field becomes 

° (a+e- mt + a-e +mt ) 

'o 



Ac 



2tt 



(5) 



Here (and throughout this paper) we approximate too ± 
f2 ~ u> inside the square root, since fl/uj ~ 3 10 -13 is 
so small; and we formally extend the integrals over Q to 
infinity, for ease of notation. 

Because the radiation pressure in the optical cavi- 
ties produces squeezing, and because this ponderomotive 
squeezing is central to the operation of our interferom- 
eters, we shall find it convenient to think about the in- 
terferometer not in terms of the single-photon modes, 
whose annihilation operators are a+ and a_, but rather 
in terms of the correlated two-photon modes (Appendix 
A and Refs. p6|,p^|) whose field amplitudes are 



a i 



a + + a_ 



V2 



a 2 



V2i 



(6) 



The commutation relations (Q) imply the following values 
for the commutators of these field amplitudes and their 
adjoints: 



[ai,a\, 



[a2,a\, 



i2w5(n - fi') 



(7a) 



and all others vanish (though some would be of order 
(fl/uj ) if we had not approximated oj ± ~ lu inside 
the square root in Eq. (5); cf. |p6|,p7|): 



[ai,oi'] = [ai,a{,] = [ 



a[,a\, 



[ai,a 2 >] 



[a], at] = , 



(7b) 



and similarly with 1 2. In terms of these two-photon 
amplitudes, Eq. (|) and = imply that the full 

electric field operator for the dark-port input is 
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TABLE I. Interferometer parameters and their fiducial val- 



Parameter 



Symbol Fiducial Value 



1.8 x 10 1J s" 



2vr x 100s~ 

30 kg 
4 km 



light frequency lu 

arm cavity ^-bandwidth 7 

grav'l wave frequency Q 

mirror mass m 

arm length L 

light power to beam splitter I 

light power to reach SQL Isql 1.0 x 10 4 W 

grav'l wave SQL /i SQL 2 x 10- 24 (-y/Q)Uz- 1/2 

opto-mech'l coupling const K, ^rfspp+lffj 

fractional signal-power loss e* 0.01 



max power squeeze factor 



0.1 



(+) 



E, 



(-) 



4n~huj 



cos(uj t) / \a\e 1 



at + a\e +iQt 



2^ 



sin(a; £) 



a 2 e + a 2 e 



dfl 
2^ 



(8) 



Thus, we see that a\ is the field amplitude for photons in 
the cos u> t quadrature and a 2 is that for photons in the 
sinc<j i quadrature [p6| , p7| . These and other quadratures 
will be central to our analysis. 

The output field at the beam splitter's dark port is de- 
scribed by the same equations as the input field, but with 
the annihilation operators a replaced by b; for example, 



4n~huj 
Ac 



cos(uj t) r (b ie - int + b\e +mt ^ 



dQ 
2^ 



+ sin(w t) / (b 2 e- lUt + b\e +int ) 



dn 
"2T 



(9) 



We shall find it convenient to introduce explicitly the 
cosine and sine quadratures of the output field, E\(t) 
and E2(i), defined by 



£out = Ex(t) cos(w Q £) + E 2 (t) s'm(uj t) 



Ej(t) = 



Ac 



dn 

2^ 



(10) 



B. Interferometer arms and gravitational waves 

LIGO's interferometers are generally optimized for 
the waves from inspiraling neutron-star and black-hole 
binaries — sources that emit roughly equal power into all 
logarthmic frequency intervals AO/fi ~ 1 in the LIGO 
band - 10Hz < / = S1/2tt < 1000Hz. Optimization 
turns out to entail making the lowest point in the inter- 
ferometer's dimensionless noise spectrum / x Sh(f) as 



low as possible. Because of the relative contributions of 
shot noise, radiation pressure noise, and thermal noise, 
this lowest point turns out to be at / = Q/2tt ~ 100 
Hz. To minimize the noise at this frequency, one makes 
the end mirrors of the interferometer's arm cavities (Fig. 
P as highly reflecting as possible (we shall idealize them 
as perfectly reflecting until Sec. VI), and one gives their 
corner mirrors transmisivities T ~ 0.033, so the cavities' 
half bandwidths are 
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Tc 
4L 



2tt x 100Hz 



(11) 



Here L = 4 km is the cavities' length (the interferometer 
"arm length"). We shall refer to 7 as the interferometer's 
optimal frequency, and when analyzing QND interferom- 
eters, we shall adjust their parameters so as to beat the 
SQL by the maximum possible amount at O = 7. In 
Table [J we list 7, L and other parameters that appear 
extensively in this paper, along with their fiducial nu- 
merical values. 

In this and the next few sections we assume, for sim- 
plicity, that the mirrors and beam splitter are lossless; 
we shall study the effects of losses in Sec. [vi| below. We 
assume that the carrier light (frequency uj ) exites the 
arm cavities precisely on resonance. 

We presume that all four mirrors ( "test masses" ) have 
masses m ~ 30 kg, as is planned for LIGO-II. 

We label the two arms n for north and e for east, and 
denote by X n and X e the changes in the lengths of the 
cavities induced by the test-mass motions. We denote by 



X = X n - X„ 



(12) 



the changes in the arm-length difference, and we regard 
a; as a quantum mechanical observable (though it could 
equally well be treated as classical Q). In the absence 
of external forces, we idealize x as behaving like a free 
mass (no pendular restoring forces). This idealization 
could easily be relaxed, but because all signals below ~ 
10 Hz are removed in the data analysis, the pendular 
forces have no influence on the interferometer's ultimate 
performance. 

The arm-length difference evolves in response to the 
gravitational wave and to the back-action influence of 
the light's fluctuating radiation pressure. Accordingly, 
we can write it as 



x(t) 



Po 

to/4 



t + 



/+00 
(Lh + x B A)e- mt 
-OO 



dn 
"2T 



(13) 



Here x is the initial value of x when a particular segment 
of data begins to be collected, p a is the corresponding 
initial generalized momentum, to/4 is the reduced mass^j 
associated with the test-mass degree of freedom x, h is 
the Fourier transform of the gravitational- wave field 



In each arm of the interferometer, the quantity measured 
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h(t) = 



+00 



iQt 



dn 
2^ 



(14) 



and Xba is the influence of the radiation-pressure back 
action. Notice our notation: x, iba and h are the in- 
dependent Fourier transforms of x(t), ccba(^) and h(t). 

Elsewhere Q] we discuss the fact that x and p a influ- 
ence the interferometer output only near zero frequency 
f2 ~ 0, and their influence is thus removed when the out- 
put data are filtered. For this reason, we ignore them 
and rewrite x{t) as 



x(t) = 



(Lh- 



XBA)e — 
Zn 



(15) 



C. Output field expressed in terms of input 

Because we have idealized the beam splitter as in- 
finitesimally thin, the input field emerging from it and 
traveling toward the arm cavities has the coherent laser 
light in the same cos uj Q t quadrature as the dark-port field 
amplitude a±. We further idealize the distances between 
the beam splitter and the arm-cavity input mirrors as 
integral multiples of the carrier wavelength A = 2irc/u) 
and as small compared to 2ttc/j ~ 300m. (These ide- 
alizations could easily be relaxed without change in the 
ultimate results.) 

Relying on these idealizations, we show in Appendix [b| 
that the annihilation operators bj for the beam splitter's 
output quadrature fields Ej(t) are related to the input 
annihilation operators a,j and the gravitational-wave sig- 
nal h by the linear relations 



bi = Abi = die 



2iii 



Abo 



'2KZ- 



ISQL 



Abo 



(a 2 -/Cai)e 2l/3 . (16) 



Here and below, for any operator A, AA = A— (A). This 
input-output equation and the quantities appearing in it 
require explanation. 

The quantities Abj are the noise-producing parts of 
bj, which remain when the gravitational- wave signal is 
turned off. The ctj impinge on the arm cavities at a 
frequency ui a + ft that is off resonance, so they acquire 
the phase shift 2/3 upon emerging, where 



f3 = arctan(f2/7) 



(17) 



is the difference between the positions of the two mirrors' 
centers of mass; this degree of freedom behaves like a free 
particle with reduced mass m r = m x m/(m + m) — m/2. 
The interferometer output is the difference, between the two 
arms, of this free-particle degree of freedom; that difference 
behaves like a free particle with reduced mass m r /2 = m/4. 



If the test masses were unable to move, then Abj would 
just be a^e 2 *' 3 ; however, the fluctuating light pressure 
produces the test-mass motion xba, thereby inducing a 
phase shift in the light inside the cavity, which shows 
up in the emerging light as the term —JCa\ in b 2 . (cf. 
Appendix pf) . The quantity 



K 



(J /Jsql)2 7 4 

ft 2 ( 7 2 + n 2 ) 



(18) 



is the coupling constant by which this radiation-pressure 
back-action converts input a\ into output Ab 2 . In this 
coupling constant, /sql is the input laser power required, 
in a conventional interferometer (Sec. HI), to reach the 
standard quantum limit: 



mL 2r y , 
/sql = —r^- 1.0 x 10 4 W . 
4w G 



(19) 



In Eq. (16), the gravitational- wave signal shows up as 
the classical piece V2KZh/hsQL of 62- Here, as we shall 
see below, 



1 SQL 



(20) 



is the standard quantum limit for the square root of the 
single-sided spectral density of h(t), y/Sh- 



III. CONVENTIONAL INTERFEROMETER 

In an (idealized) conventional interferometer, the 
beam-splitter's output quadrature field is measured 
by means of conventional photodetection.Q The Fourier 
transform of this measured quadrature is proportional 
to the field amplitude b 2 = Ab 2 + v2K.(h / '/isQL)e I/3 '; cf. 
Eqs. (|Io| ) and (|l^). Correspondingly, we can think of 
62 = ^2^) as the quantity measured, and when we com- 
pute, from the output, the Fourier transform h = h(£l) 
of the gravitational-wave signal, the noise in that com- 
putation will be 



4 Here and throughout this paper we regard some particu- 
lar quadrature E((t) as being measured directly. This corre- 
sponds to superposing on E((t) carrier light with the same 
quadrature phase as E( and then performing direct photode- 
tection, which produces a photocurrent whose time variations 
are proportional to E((t). For a conventional interferometer 
the carrier light in the desired quadrature, that of i?2(i), can 
be produced by operating with the dark port biased slightly 
away from the precise dark fringe. In future research it might 
be necessary to modify the QND designs described in this 
paper so as to accommodate the modula tions that are actu- 
ally used in the detection process; see Sec. VIII and especially 
Footnote O. 
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h n (Q) = ^Alu 



'2K. 



iff 



(21) 



This noise is an operator for the Fourier transform of a 
random process, and the corresponding single-sided spec- 
tral density Sh(f) associated with this noise is given by 
the standard formula fitMETl 



\2n8(n - n')s h (f) = (m\h n (n)hUn')H 



sym • 



(22) 



Here / = Q/2n is frequency, |in) is the quantum state of 
the input light field (the field operators a\ and 02), and 
the subscript "sym" means "symmetrize the operators 
whose expectation value is being computed", i.e., replace 
h n {Q)hl(Q,') byUh n (n)hUn') + hi(tt)h n {n)YNote 
that when Eq. ( |2l| ) for h n is inserted into Eq. (f22[), the 
phase factor e _l ^cancels, i.e. it has no influence on the 
noise Sh- This allows us to replace Eq. (EHl) by 



'2JC 



(23) 



For a conventional interferometer, the dark-port input 
is in its vacuum state, which we denote by 



I in) = |0 a ) 



(24) 



For this vacuum input, the standard relations a + |0 a ) = 
a_|0 a ) = 0, together with Eqs. (§) and (0), imply [|fj27| 



(0 o |aj(4,|0 o ) S y m = ^2n5(n- tt')8j k 



(25) 



Comparing this relation with Eq. (|22|) and its generaliza- 
tion to multiple random processes, we see that (when 
|in) = |0 a )) ai(Q) and 02(0) can be regarded as the 
Fourier transforms of classical random processes with 
single-sided spectral densities and cross-spectral density 
given by [g 



S ai (f)=S a2 (/)=!, S aia2 (f) = 



(26) 



Combining Eqs. © and ©-(H) [or, equally well, 
(pf), (p3|), and (p6|)], we obtain for the noise spectral 
density of the conventional interferometer 



S h = 



h 2 

"-SQL 



(27) 



This spectral density is limited, at all frequencies Q, by 
the standard quantum limit 



> h 2 

±- "-SQL 



STi 



in 2 L 2 



(28) 



Recall that IC is a function of frequency £1 and is pro- 
portional to the input laser power I [Eq. In our 
conventional interferometer, we adjust the laser power to 
Io = /sql [Eq. ©], thereby making /C(0 = 7) = 1, 




squeezed-variational 



0.1 0.2 



0.5 



1 

QJy 



10. 



FIG. 4. The square root of the spectral density \/Sh of the 
gravitational-wave noise for several interferometer designs, as 
a function of angular frequency Q, with optical losses assumed 
negligible; \fSh is measured in units of the standard quantum 
limit at frequency Q — 7, and f2 is measured in units of 7. The 
noise curves shown are: (i) the standard quantum limit itself, 
/isql(£2) [Eq. (|20|]; (ii) the noise for a conventional interfer- 
ometer with laser power I a — Isql [Eq. (^9|)] ; (iii) the noise for 
a squeezed-input interferometer with I — Jsql, squeeze factor 
e~ 2R — 0.1, and (a) optimized FD squeeze angle A = — $(fi) 
[Eq. (|49|); solid curve] , (b) optimized frequency-independent 
squeeze angle [Eq. (|52|); dashed curve]; (iv) the noise for a 
variational- output interferometer with I — 10/sql and opti- 
mized frequency- dependent homodyne phase ( = 3>(fi) [Eq. 
(S)]; and (v) the noise for a squeezed-variational interferom- 
eter with I — 10/sql, input squeeze factor e~ 2R — 0.1, and 
optimized input squeeze angle X — tt/2 and output homodyne 
phase C = [Eq. ®]. 

which minimizes Sh at the interferometer's optimal fre- 
quency Q = 7. The noise spectral density then becomes 
[cf. Eqs. (jH) and ©] 



Sh 



Ah 



iL 2 n 2 



2 7 4 



fi 2 (7 2 



n 2 



n 2 {-y 2 + n 2 ) 



2 7 4 



(29) 



This optimized conventional noise is shown as a curve 
in Fig. |J, along with the standard quantum limit hsQL 
and the noise curves for several QND interferometers to 
be discussed below. This conventional noise curve is cur- 
rently a tentative goal for LIGO-II, when operating with- 
out signal recycling . 



IV. STRATEGIES TO BEAT THE SQL, AND 
THEIR LOSSLESS PERFORMANCE 

A. Motivation: Ponderomotive squeezing 

The interferometer's input-output relations A61 = 
aie 2l/3 , A62 = (a 2 — K.ai)e 2zl3 can be regarded as consist- 
ing of the uninteresting phase shift e 2l/3 , and a rotation in 
the {01,02} plane (i.e., {cos oj t, sinw t} plane), followed 
by a squeeze: 



bj = S+(r, (j})R ji (-9)a J e 2i0 R{-e)S(r, 0) 



(30) 
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Here R(—0) is the rotation operator and S(r, cj>) the 
squeeze operator for two-photon quantum optics; see Ap- 
pendix |A| for a very brief summary, and Refs. (2(]j27j] for 
extensive detail. The rotation angle 9, squeeze angle <f> 
and squeeze factor r are given by 

9 = arctan(/C/2), = -arccot(/C/2), r — arcshinh(/C/2). 



Note that, because the coupling constant K, depends on 
frequency £1 [Eq. (|l8|)1, the rotation angle, squeeze angle, 
and squeeze factor are frequency dependent. This fre- 
quency dependence will have major consequences for the 
QND interferometer designs discussed below. 

The rotate-and-squeeze transformation ( |30| ) for the 
two-photon amplitudes implies corresponding rotate- 
and-squeeze relations for the one-photon creation and 
annihilation operators 



b± = S\r,4>)R ] {-9)a±e ±2lfl R(-9)S{r,ct>) 



(32) 



Denote by |0 a+ ) the vacuum for the in mode at frequency 
ll> + Q, by |0 a _) that for the in mode at lj q — VL, and by 
|0 a± ) the vacuum for one or the other of these modes; 
and denote similarly the vacuua for the out modes, |0b±). 
Then |0 o± ) is the state annihilated by a± and |0& ± ) is that 
annihilated by b±. Correspondingly, the rotate-squeeze 
relation (|32|) implies that 



b±\0 b± ) = S i R' i a ± e ±2l/3 RS\0 b± ) 







(33) 



where the parameters of the squeeze and rotation oper- 
ators are those given in Eqs. (|3l]) and (|32"|). This equa- 
tion implies that e ±2t P RS\0b ± ) is annihilated by a± and 
therefore is the in vacuum |0 a± ) for the in mode uj ± tt: 



e ±2 ^RS\0 b± ) = |0 a± ) . 

Applying i?t and noting that R^\0 a± ) 
uum is rotation invariant), we obtain 



|0» ± ) 



|06±) 



(34) 

|0 Q± ) (the vac- 



(35) 



Thus, the in vacuum is equal to a squeezed out vacuum, 
aside from an uninteresting, frequency-dependent phase 
shift. The meaning of this statement in the context of a 
conventional interferometer is the following: 

For a conventional interferometer, the in state is 



) = |O a± )=e ±2 ^(r,0)|O b± > 



(36) 



and because we are using the Heisenberg Picture where 
the state does not evolve, the light emerges from the in- 
terferometer in this state. However, in passing through 
the interferometer, the light's quadrature amplitudes 
evolve from aj to bj. Correspondingly, at the output we 
should discuss the properties of the unchanged state in 
terms of a basis built from the out vacuum \0 b± ). Equa- 
tion (|35|) says that in this owi language, the light has been 



squeezed at the angle <f> and squeeze-factor r given by Eq. 
( pl| ) . This squeezing is produced by the back-action force 
of fluctuating radiation pressure on the test masses. That 
back action has the character of a ponderomotive nonlin- 
earity first recognized by Braginsky and Manukin ]l5[|.p| 
The correlations inherent in this squeezing form the foun- 
dation for the QND interferometers discussed below. 
One can also deduce this ponderomotive squeezing 



(31) from the in-out relations A£>i = aie 2l/3 , Ab 2 = (a 2 



/Cai)e 2j/3 [Eq. (fig)], the expressions 



-2ir6{n - n')S bj (f) = (in\AbjAbjS\m) syia , 



-2nS(n-n')s blb2 (f) = 



(in|i(A6iA6 2 , t + A6 t 1 A6 2 0|in 



sym 

(37) 



for the spectral densities and cross spectral densities of 
b\ and 62, and the spectral densities S ai = S a2 = 1, 
S ai a 2 — [Eqs. (^6|)]. These imply that for a conven- 
tional interferometer 



S bl = 1 , S b2 = 1 + K? , S blb2 = —K 



(38) 



Rotating Abj through the angle cf> = iarccot(/C/2) to 
obtain 

b\ = bi cos <f) + 62 sin <f> , b' 2 = b 2 cos <j) — b\ sin (p , (39) 

and using Eqs. (|37]) and ([58]), we obtain 

2 



Sh> = e 



-2r 



Si, 



+2r 



( y/1 + (/C/2) 2 - AC/2) ~1//Cif/C>l, 
(v/l + (/C/2)2 + /C/2) 2 , S ViK =0, 



(40) 



which represents a squeezing of the input vacuum noise 
in the manner described formally by Eqs. ( |36| ) and (|3l|). 

This ponderomotive squeezing is depicted by the noise 
ellipse of Fig. || For a conventional interferometer (b 2 
measured via photodetection^j) , the signal is the arrow 
along the b 2 axis, and the square root of the noise spectral 
density S b2 is the projection of the noise ellipse onto the 
b 2 axis. For a detailed discussion of this type of graphical 
representation of noise in two-photon quantum optics see, 
e.g., Ref. M . 



5 Recently it has been recognized that this ponderomotive 
nonlinearity acting on a movable mirror in a Fabry-Perot res- 
onator may provide a practical method for generating bright 
squeezed light [3C 

6 See footnote 
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FIG. 5. Noise ellipses for a conventional interferometer. 
Left: Noise for vacuum that enters the interferometer's dark 
port. Right: Noise for ponderomotively squeezed vacuum 
that exits at the dark port along with the gravitational-wave 
signal; the ponderomotive squeeze has moved the point P 
to the new point P' [b± = ai, 62 = 02 — /Can,, Eqs. (|l6|)]. 
These noise ellipses have dimensions and shapes described by 
the noise spectral densities (H), (^) and (icj), and by the 
squeeze equations ( pioj ) a nd (|31| ). The minor radius of the 
output noise ellipse is y/B^i = e~ r , and its major radius 
is yJS b i = e +r , where r is the squeeze factor; cf. Eqs. ( j3l| ) 
and @. The conventional interferometer measures i>2 , which 
contains the indicated noise [cf. Eq . (B3 )] and the indicated 
signal [Sb-2 = V^K-h/hsQL', cf. Eq. (|16|)j. For a detailed dis- 
cussion of noise ellipses in 2- photon quantum optics see, e.g., 
Ref. H. 



B. Squeezed-Input Interferometer 

Interferometer designs that can beat the SQL (E8I) are 
sometimes called "QND interferometers" . Unruh pi has 
devised a QND interferometer design based on (i) putting 
the input electromagnetic fluctuations at the dark port 
(ai and 02) into a squeezed state, and (ii) using standard 
photodetection to measure the interferometer's output 
field. We shall call this a squeezed-input interferome- 
ter. The squeezing of the input has been envisioned as 
achieved using nonlinear crystals pO 21 1, but one might 
also use ponderomotive squeezing. 

The squeezed-input interferomete r is identical to the 
conventional interferometer of Sec. 



Ill, except for the 



choice of the in state |in) for the dark-port field. Whereas 
a conventional interferometer has |in) = |0 o ), the 
squeezed-input interferometer has 



|in) =S(R,X)\O a ) 



(41) 



where R is the largest squeeze factor that the experi- 
menters are able to achieve (e~ 2R ~ 0.1 in the LIGO-III 
time frame), and A = X(Q) is a squeeze angle that de- 
pends on side-band frequency. One adjusts X(Q) so as to 
minimize the noise in the output quadrature amplitude 
62, which (i) contains the gravitational- wave signal and 
(ii) is measured by standard photodetection. 

The gravitational- wave noise for such an interferometer 
is proportional to 



(0 a \h r , 



|o«) 



(42) 



[Eq. (|22])], where h ns is the squeezed gravitational-wave 
noise operator 



Ks = S*(R,X)h n S(R, X) 



(43) 



and h' n = h n (Q'). By inserting expression (|2l| ) for h n 
into Eq. ( [43] ) and then combining the interferometer's 
ponderomotive squeeze relation A62 = (02 — /Cai)e 2l/3 
with the action of the squeeze operator on a\ and 02 [Eq. 
|)], we obtain 



*SQL 

/2X 



v/(l + /C 2 )e 4/3 



x (at {cosh R cos $ - sinhi?cos[$ - 2(<I> + A)]} 
- a 2 { cosh i? sin $ - sinhi?sin[$ - 2($ + A)]}) , (44) 



where 



$ = arccot/C . 



(45) 



We can read the spectral density of the gravitational- 
wave noise off of Eq. (^) by recalling that in the |0 a ) 
vacuum state [which is relevant because of Eq. p^)], a\ 
and ai can be regarded as random processes with spec- 
tral sensities S a -, — S a2 — 1 and vanishing cross spectral 
density [Eqs. (fa)]: 



Si, 



h 2 

SQL 



K) (cosh2i?-cos[2(A + $)]sinh2i? 

(46) 



It is straightforward to verify that this noise is minimized 
by making it proportional to cosh2i? — sinh2i? = e~ 2R , 
which is achieved by choosing for the input squeeze angle 



A(fi) = -$(0) = -arccot/C(fi) 
The result is 



Sh - — {lc +IC 



-2R 



(47) 



(48) 



This says that the squeezed-input interferometer has 
the same noise spectral density as the conventional inter- 
ferometer, except for an overall reduction by e~ 2R , where 
R is the squeeze factor for the dark-port input field (a re- 
sult deduced by Unruh |lj] and later confirmed by Jaekel 
and Reynaud jl8| using a different method); see Fig. |^. 
This result implies that the squeezed-input interferome- 
ter can beat the amplitude SQL by a factor e~ R . 

When the laser power I a of the squeezed-input interfer- 
ometer is optimized for detection at the frequency fl = 7 
(I = Isql a s for a conventional interferometer), the 
noise spectrum becomes 



Sh 



Ah 



iL 2 n 2 



27 4 



n 2 (-f 2 + n 2 ) 



Q 2 {j 2 + Q 2 ) 



2 7 4 



-2ft 



(49) 
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This optimized noise is shown in Fig. || for e~ 2R = 0.1, 
along with the noise spectra for other optimized interfer- 
ometer designs. 

In previous discussions of this squeezed-input scheme 
| p^ , ^8|JT9| , no attention has been paid to the practical 
problem of how to produce the necessary frequency de- 
pendence 



A(O) = = -arccot 



27 4 



ft 2 ( 7 2 + Q 2 



(50) 



of the squeeze angle. In Sec. |V Q , we shall show that 
this A(f2) can be achieved by squeezing at a frequency- 
independent squeeze angle (using, e.g., a nonlinear 
crystal for which the squeeze angle will be essentially 
frequency-independent because the gravity-wave band- 
width, < 1000 Hz, is so small compared to usual optical 
bandwidths ) and then filtering through two Fabry-Perot 
cavities. This squeezing and filtering must be performed 
before injection into the interferometer's dark port; see 
Fig. [l] for a schematic diagram. 

The signal and noise for this squeezed-input interfer- 
ometer are depicted in Fig. ^|. 

We comment, in passing, on two other variants of a 
squeezed-input interferometer: 

(i) If, for some reason, the filter cavities cannot be 
implemented successfully, one can still inject squeezed 
vacuum at the dark port with a frequency-independent 
phase that is optimized for the lowest point in the noise 
curve, f2 = 7; i.e. (with the input power optimized to 
lo = Isql)- 



A = -$(7) = — 7r/4 ; 



(51) 



cf. Eq. (p0|). In this case the noise spectrum is 
S , = %(i + K 



[(cosh R cos $ - sinh R sin $) 2 
h(coshi?sin$ - sinh R cos $) 2 ] 



(52) 



^ = 0.1. The SQL is beat by the same 
2R ~ 0.32 as in the case of a fully opti- 



[Eq. ([46]), translated into gravitational- wave noise via Eq. 
(p3)]. This noise spectrum is shown as a dashed curve 
in Fig. ||, for e 
factor fx — \[e 
mized squcczed-input interferometer, but the frequency 
band over which the SQL is beat is significantly smaller 
than in the optimized case, and the noise is worse than 
for a conventional interferometer outside that band. 

(ii) Caves in a paper that preceeded Unruh's 

and formed a foundation for Unruh's ideas, proposed 
a squeezed-input interferometer with the squeeze angle 
set to A = 7r/2 independent of frequency. In this case, 



Eq. (46), translated into gravitational- wave noise via Eq. 
(p3|), says that 



Sh 



h 2 

"•SQL 



3 2 «/C 



, 2R X 



(53) 



a 2' 




• 
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FIG. 6. Noise ellipses for a squeezed-input interferometer. 
Left: Noise for squeezed vacuum that enters the interferom- 
eter's dark port. The field is squeezed at the angle A = — <fr. 
Right: Noise for the field that exits at the dark port along 
with the gravitational-wave signal. This output field results 
from the interferometer's ponderomotive squeezing of the in- 
put field [e.g., point P goes to point P' in accord with 61 =01, 
&2 = C12 — /Can,; Eqs. (JTsh] . If the input field had been vacuum 
as in a conventional interferometer [Fig. [sj, then the output 
would have been squeezed in the manner of the dashed el- 
lipse. The two squeezes (input and ponderomotive) result in 
the shaded ellipse, whose projection along the axis measured 
by the photodetector (62 axis) has been minimized by the 
choice of squeeze angle, A = — $. 



Since K, is proportional to the input laser power I a , Caves' 
interferometer produces the same noise spectral density 
as a conventional interferometer [Eq. (|27|)] but with an 
input power that is reduced by a factor e~ 2R . This is a 
well-known result. 



C. Variational-output interferometer 

Vyatchanin, Matsko and Zubova p^-|25| have devised a 
QND interferometer design based 011(1) leaving the dark- 
port input field in its vacuum state, |in) = |0 Q ), and (ii) 
changing the output measurement from standard pho- 
todetection (measurement of 62) to homodyne detection 
at an appropriate, frequency-dependent (FD) homodyne 
phase C(f2) ~~ i- e -? measurement of 



= 61 cos £ + 62 sin ( . 



(54) 



In their explorations of this idea, Vyatchanin, Matsko 
and Zubova ]23|-p5[| did not identify any practical scheme 
for achieving such a FD homodyne measurement, so they 
approximated it by homodyne detection with a homo- 
dyne phase that depends on time rather than frequency 
— a technique that they call a "quantum variational mea- 
surement" . 

In Sec. below, we show that the optimized FD homo- 
dyne measurement can, in fact, be achieved by filtering 
the interferometer output through two Fabry-Perot cav- 
ities and then performing standard, balanced homodyne 
detection at a frequency-independent homodyne phase; 
see Fig. || for a schematic diagram. We shall call such an 
scheme a variational-output interferometer. The word 
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"variational" refers to (i) the fact that the measurement 
entails monitoring a frequency- varying quadrature of the 
output field, as well as (ii) the fact that the goal is to mea- 
sure variations of the classical force acting on the inter- 
ferometer's test mass (the original Vyatchanin-Matsko- 
Zubova motivation for the word). 

The monitored FD amplitude [Eq. (pi])] can be ex- 
pressed in terms of the interferometer's dark-port in- 
put amplitudes a±, a 2 and the Fourier transform of the 
gravitational- wave field h as 



sinC V2/C 



n SQL 



+ [a2 + (cot( - )C)ai]e 



2ift 



(55) 

cf. Eqs. Jl6| ) and (|54|). Correspondingly the operator 
describing the Fourier transform of the interferometer's 
gravitational-wave noise is 



hn(Q) = -p^e l/3 [a 2 + ai(cotC-/C)] 



2AC 



(56) 



cf. Eq. (f23|). 

The radiation-pressure-induced back action of the 
measurement on the interferometer's test masses is em- 
bodied in the —Ka\ term of this equation; cf. Eq. Jl6| ) 
and subsequent discussion. It should be evident that by 
choosing 



C = <f> = arccot/C 



(57) 



we can completely remove the back-action noise from the 
measured interferometer output; cf. Fig. 0. This optimal 
choice of the FD homodyne phase, together with the fact 
that the input state is vacuum, |in) = |0 a ), leads to the 
gravitational- wave noise 



S h = 



h 2 1 
SQL 1 



1 / Ah \n 2 (j 2 + n 2 ) 



2 K Io/Isql \mL 2 n 2 



2 7 4 



(58) 



cf. Eqs. (|22j) and @. 

This noise for an optimized variational-output inter- 
ferometer is entirely due to shot noise of the measured 
light, and continues to improve cx 1/I even when the 
input light power I a exceeds Isql- Figure f| shows this 
noise, along with the noise spectra for other optimized 
interferometer designs. 

It is interesting that the optimal frequency-dependent 
homodyne phase $ for this variational-output inter- 
ferometer is the same, aside from sign, as the opti- 
mal frequency-dependent squeeze angle for the squeezed- 
input interferometer; cf. Eq. (E^). 



D. Comparison of squeezed-input and 
variational-output interferometers 

The squeezed-input and variational-output interferom- 
eters described above are rather idealized, most especially 
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FIG. 7. Noise ellipses for a variational-output interferome- 
ter. Left: Noise for the ordinary vacuum that enters the inter- 
ferometer's dark port. Right: Noise for the field that exits at 
the dark port along with the gravitational-wave signal. These 
noise ellipses are the same as for a conventional interferome- 
ter, Fig. g, but here the quantity measured is the quadrature 
amplitude 6$ with frequency dependent phase $ = arccot/C. 
It is informative to compare the measured phase $ with the 
angle of ponderomotive squeeze <j> = iarccot(/C/2). They are 
related by tan $ = | tan 2<j> = tan </>/(l — tan 2 <f>), so $ is al- 
ways larger than <j>; but for large K, (strong beating of the 
SQL), they become small and nearly equal. 



because they assume perfect, lossless optics. When wc re- 
lax that assumption in Sec. VI below 
for realistic squeeze factors e 



we shall see that, 
R and losses e*, the two 



interferometers have essentially the same performance, 
but the variational-output intefcrmomctcr requires ~ 10 
times higher input power I Q . In this section we shall 
seek insight into the physics of these interferometers by 
comparing them in the idealized, lossless limit. 

Various comparisons are possible. The noise curves 
in Fig. H illustrate one comparison: When the FD ho- 
modyne angle has been optimized, a lossless variational- 
output interferometer reduces shot noise below the SQL 
and completely removes back-action noise; by contrast, 
when the FD squeeze angle has been optimized, a 
squeezed-input interferometer reduces shot noise and re- 
duces but does not remove back-action noise; cf. Eqs. 
©and®. 

In variational-output interferometers, after optimizing 
the FD homodyne angle, the experimenter has further 
control of just one input/output parameter: the laser 
intensity or equivalently 7 /Jsql = JC(ft = 7). When 
Io/Isql is increased, the shot noise decreases; inde- 
pendent of its value, the back-action noise has already 
been removed completely; cf. Eq. (pq). By contrast, in 
squeezed-input interferometers, after optimizing the FD 
squeeze phase, the experimenter has control of two pa- 
rameters: Io/Isql, which moves the minimum of the 
noise curve back and forth in frequency but does not 
lower its minimum ]17j] , and the squeeze factor R, which 
reduces the noise by e~ 2R ; cf. Eq. (48). 

Present technology requires that R be approximately 
constant over the LIGO frequency band. However, in the 
same spirit as our assumption that the FD homodyne 
phase can be optimized at all frequencies, it is instruc- 
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tive to ask what can be achieved with an unconstrained, 
frequency-dependent (FD) squeeze factor R(fl), when 
coupled to an unconstrained FD squeeze angle A(f2). 

One instructive choice is A = —arccot/C as in our pre- 
vious, optimized interferometer [Eq. (ft7|)], and e~ 2R — 
1/(1 + /C 2 ). In this case, the squeezed-input interferome- 
ter has precisely the same noise spectrum as the lossless 
variational-output interferometer 



Sh 



h 2 

"-SQL 

2/C 



(59) 



[Eq. (pq)], and achieves it with precisely the same laser 
power. 

Another instructive choice is an input squeeze that is 
inverse to the interferometer's ponderomotive squeeze (a 
configuration we shall call "inversely input squeezed" or 
IIS): Let the dark-port input field before squeezing be 
described by annihilation operators c±, so 



c±|m) 



Q 



(60) 



i.e. the pre-squeeze field is vacuum. Then, denoting by 
ci , c 2 the quadrature amplitudes of this pre-squeeze field, 
the IIS input squeezing is 



ai = Ci 



«2 



c 2 + /Cci , 



(61) 



where /C(f2) is the interferometer's frequency-dependent 
coupling constant (p^). The interferometer's ponderomo- 
tively squeezed output noise is then 

Ah = a x e 2lfi = c x e 2lf> , Ab 2 = (a 2 - JCa x )e 2lP = c x e 2lf3 

(62) 

[cf. Eq. (|l6|)1, i.e., the noise of the output light is that 
of the vacuum with a phase shift, but since the vacuum 
state is insensitive to phase, it is actually just the noise 
of the vacuum. 

Figure H illustrates this: The IIS input light is squeezed 
in a manner that gets perfectly undone by the pondero- 
motive squeeze, so the output light has no squeeze at all. 
The fact that the input squeeze is inverse to the pon- 
deromotive squeeze shows up in this diagram as an input 
noise ellipse that is the same as the output ellipse of the 
ponderomotively squeezed vacuum, Fig. ||, except for a 
reflection in the horizontal axis. 

Because the output of the IIS interferometer is b 2 (or- 
dinary photodetection) and the output light's state is the 
ordinary vacuum, its gravitational-wave noise is 
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2/C 



(63) 



cf. Eqs. (||), (||) and (§|) (with a, replaced by b 3 ). No- 
tice that this is identically the same noise spectral den- 
sity as for our previous example [Eq. J59|)] and as for a 
variational-output interferometer, and it is achieved in 
all three cases with the same light power. 
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FIG. 8. Noise ellipses for a squeezed-input interferometer 
whose input squeeze is inverse to the interferometer's pon- 
deromotive squeeze ("IIS interferomter" ) . 



The fact that our two squeezed-input examples pro- 
duce the same noise spectrum using different squeeze an- 
gles and squeeze factors should not be surprising. The 
noise spectrum is a single function of f2 and it is being 
shaped jointly by the two squeeze functions A(f2) and 

The fact that the IIS interferometer and the varia- 
tional output interferometer produce the same noise spec- 
tra results from a reciprocity between the IIS and the 
variational-output configurations: The IIS interferometer 
has its input squeezed at the angle — $ = —arccot/C and it 
has vacuum-noise output, whereas the variational-output 
interferometer has vacuum-noise input and is measured 
at the homodyne angle +$ = +arccot/C. 

Note that the IIS interferometer has a different input 
squeeze ang le [A = -±arccot(/C/2); cf. Eq. @] from 
that of the angle- optimized squeezed- input interferome- 
ter of Sec. [iVBl [A = -arccot/C; cf. Eq. @]. This dif- 
ference shows clearly in the noise ellipses of Fig. || (the 
IIS interferometer) and Fig. ^ (the angle-optimized in- 
terferomter). Meoreover, this difference implies that by 
optimizing the IIS interferometer's squeeze angle (chang- 
ing it to A = —arccot/C), while keeping its squeeze factor 
unchanged [R — arcshinh(/C/2); cf. Eq. (|3l|)1, we can im- 
prove its noise performanance slightly. The improvement 
is from (63) to 



S h = 
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1 + /C 2 



1 + \ (JC 2 + /CV/C 2 + 4) 



(64) 



[which can be derived by setting A = — $ = —arccot/C 
and R = arcshinh(/C/2) in Eq. ([Iq), or by inserting R = 
arcshinh(/C/2) into Eq. @— note that @ is valid for 
any angle-optimized, squeezed-input interferometer but 
not for the IIS interferometer]. The improvement factor 
in square brackets is quite modest; it lies between 0.889 
and unity. 

We reiterate, however, that the above comparison of 
interferometer designs is of pedagogical interest only. In 
the real world, the noise of a QND interferometer is 
strongly influenced by losses, which we consider in Sec. 
VI below. 
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E. Squeezed-variational interferometer 

The squeezed-input and variational-output techniques 
are complementary. By combining them, one can beat 
the SQL more strongly than using either one alone. We 
call an interferometer that uses the two techiques simul- 
taneously a squeezed-variational interferometer. 

The dark-port input of such an interferometer is 
squeezed by the maximum achievable squeeze factor R 
at a (possibly frequency dependent) squeeze angle A(fi), 
so 



S(R,X)\O a ) 



(65) 



The dark-port output is subjected to FD homodyne de- 
tection with (possibly frequency dependent) homodyne 
angle C(^)i i- e -j the measured quantity is the same out- 
put quadrature as for a variational-output interferome- 
ter, [Eq. (p5|)], so the gravitational- wave noise operator 
is also the same 

h n (fl) = e lli [oa + ai(cotC- tC)]e lp (66) 



/2K, 
[Eq. ©]. 

As for 
a squeezed-input interferometer, the gravitational-wave 
noise is proportional to 



(in | /i„/i n / 1 in) = {O a \h ns h ns > |0 Q ) 



(67) 



[Eq. (|22|)], where h ns is the squeezed gravitational- wave 
noise operator 



S*{R, X)h n S(R, A) 



(68) 



By inserting expression (|66J) for h n into Eq. (|68|) and 
invok ing the action of the squeeze operator on a\ and 
[Eq. flAq)], we obtain 
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x ^ai{coshi?cosl> - sinhi?cos[«l - 2(<l + A)]} 
-a2{coshi?sin«i - sinhi?sin[l> - 2(l> + A)]} 
where 

JC = JC — cot £ , $ = arccotAC . 



(69) 



(70) 



As for a squeezed-input interferometer [see passage fol- 
lowing Eq. (Hf|)], we can read the gravitational- wave spec- 
tral density off of Eq. (|69| ) by regarding a\ and a-i as ran- 
dom processes with unit spectral densities and vanishing 
cross spectral density. The result is 
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FIG. 9. Noise ellipses for a squeezed-variational interfer- 
ometer. Left: Noise for the squeezed vacuum that enters the 
interferometer's dark port. Right: Noise for the field that ex- 
its at the dark port along with the gravitational- wave signal. 



This noise is minimized by setting the input squeeze 
angle A and output homodyne phase £ to 



A = 7r/2 , ( = $ = arccot/C , 
which produces JC — and A = $ = 7r/2, so 
h 2 

a SQL -2R 
oh = — e 



(72) 
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(73) 



V^sql \mL 2 SP, 
see Fig. | 

Equation ( |7^ ) says that, to optimize the (lossless) 
squeezed-variational interferometer, one should squeeze 
the dark-port input field at the frequency-independent 
squeeze angle C = 77 /2 (which ends up squeezing the 
interferometer's shot noise), and measure the output 
field at the same FD homodyne phase ( — $ as for a 
variational-output interferometer; see Fig. ^. Doing so 
produces an output, Eq. (|73]), in which the radiation- 
pressure-induced back-action noise has been completely 
removed, and the shot noise has been reduced by the 
input squeeze factor e~ 2R . 

Because the optimal input squeeze angle is frequency 
independent, the squeezed variational interferometer 
needs no filter cavities on the input. However, they are 
needed on the output to enable FD homodyne detection; 
see Fig. || for a schematic diagram. 



V. FD HOMODYNE DETECTION AND 
SQUEEZING 

Each of the QND schemes discussed above requires 
homodyne detection with a frequency-dependent phase 
(FD homodyne detection) and/or input squeezed vacuum 
with a frequency-dependent squeeze angle (FD squeezed 
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vacuum). In this section we sketch how such FD homo- 
dyne detection and squeezing can be achieved. 



If we adjust the filter and the constant homodyne 
phase so that 



A. General Method for FD Homodyne Detection 

The goal of FD Homodyne Detection is to measure the 
electric-field quadrature 



E c (t) 



(V 



-iQt 



Ac jq 

for which the quadrature amplitude is 

= bi cos £ + 6 2 sin C , C — C(^) 



bU +iUt )^, (74) 
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cf. Eqs. ( p~0| ) and (H). If ( were frequency indepen- 
dent, the measurement could be made by conventional 
balanced homodyne detection, with homodyne phase £. 
In this subsection we shall show that, when £ depends on 
frequency, the measurement can be achieved in two steps: 
first send the light through an appropriate filter ( assumed 
to be lossless), and then perform conventional balanced 
homodyne detection. 

The filter puts onto the light a phase shift a that de- 
pends on frequency. Let the phase shift be a+ for light 
frequency u> + f2, and a_ for ujq — tt. The input to the 
filter has amplitudes (annihilation operators) b± at these 
two sidebands, and the filter output has amplitudes (de- 
noted by a tilde) 

b± = b±e m± . (76) 
The corresponding quadrature amplitudes are 



h = 



V2 



h = 



b+ - bl 
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(77) 



at the input [Eqs. @], and the analogous expression with 
tildes at the output. Combining Eqs. ( |77| ) with and 
without tildes, and Eq. (|76|), we obtain for the output 
quadrature amplitudes in terms of the input 



bi = e IQm (&i cosa p — 62 sina p ) , 
&2 = e JQm (62 cos a p + 62 sin a p ) . 



Here 



-{a+ +a_) 



(78) 



(79) 



The light with the output amplitudes bi , 62 is then sub- 
jected to conventional balanced homodyne detection with 
frequency-independent homodyne angle 9, which mea- 
sures an electric-field quadrature with amplitude 



bg = bi cos + 62 sin 9 

= e la,T1 [bi cos(# - a p ) + b 2 sin( 



(80) 



-(a + +a_) = C(0) , 



(81) 



then, aside from the frequency-dependent phase shift a m , 
the output quadrature amplitude will be equal to our 
desired FD amplitude: 



b e = 



(82) 



The phase shift a m (Q) is actually unimportant; it can 
be removed from the signal in the data analysis (as can 
be the phase shift produced by the interferometer's 
arm cavities). 

To recapituate: FD homodyne detection with homo- 
dyne phase C(^) can be achieved by filtering and conven- 
tional homodyne detection, with the filter's phase shifts 
a± (at uj = uj ±H) and the constant homodyne phase 9 
adjusted to satisfy Eq. 



B. Realization of the filter 



The desired FD homodyne phase is 
( = ^(fi) = arccot/C = arccot 
= arctan 



A 4 



n 2 



A 4 



where 



A 
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(83) 



(84) 



[cf. Eqs. flU) and @]. Recall that 7 ~ 2ir x 100Hz is 
the optimal frequency of operation of the interferometer, 
and to beat the SQL by a moderate amount will require 
Io/Isql - 10 so A 4 ~ 20 7 4 , i.e., A ~ 2 7 . 

In Appendix ^| we show that this desired FD phase can 
be achieved by filtering the light with two successive loss- 
less Fabry-Perot filter cavities, followed by conventional 
homodyne detection at homodyne angle 



= tt/2 



(85) 



[i.e., homodyne measurement of 62 at the filter output; 
cf. Eq. (§(])].[] The two filter cavities (denoted I and II) 



7 The fact that only two cavities are needed to produce the 
desired FD homodyne phase (^) is a result of the simple 
quadratic form of tan<E>(f2 2 ). If the desired phase were sig- 
nificantly more complicated , a l arger number of filter cavities 
would be needed; cf. Eq. (S3) and associated analysis. It 
would be interesting to explore what range of FD homodyne 
phases can be achieved, with what accuracy, using what num- 
ber of cavities. 
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produce phase shifts ctj± and an± on the ui ± f2 side 
bands, so upon emerging from the second cavity, the net 
phase shifts are 



a± = a i± + a u ± 



(86) 



Each cavity (J =1 or II) is characterized by two pa- 
rameters: its decay rate (bandwidth) 28j (with J 
1 or II), and its fractional resonant-frequency offset from 
the light's carrier frequency io , 



Wres J 



'J 



(87) 



Here w res j is the resonant frequency of cavity J. In terms 
of these parameters, the phase shifts produced in the 
lo ± 51 side bands by cavity J are 



cn,/± = arctan(£j ± Cl/Sj) 



(88) 



The filters' parameters must be adjusted so that the net 
phase shift (|S6"|), together with the final homodync angle 
9 = 7r/2, produce the desired FD phase, Eqs. (|8l| ) and 

In Appendix ^ we derive the following values for the 
filter parameters £i, Si, £n and Su as functions of the 
parameters A and 7 that appear in the desired FD ho- 
modyne phase. Define the following four functions of A 
and 7 



A + = 



Q + VQ 
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Then in terms of these functions, the filter parameters 
are 



(89c) 
(89d) 
(89e) 
(89f) 
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Note that, when the cavity half-bandwidths 5j are 
expressed in terms of the half-bandwid th 7 of th e in - 
terferometer's arm cavities, as in Eqs. (89e) and (891), 
then the filter parameters depend on only one character- 
istic of the desired FD homodyne phase: the quantity 
(A/7) 4 = 2/ //sql- Figure [l0| depicts the filter parame- 
ters as functions of this quantity. 

As Fig. [l(] shows, the half-bandwidths of the two filter 
cavities are within a factor ~ 2 of that of the interferom- 
eter's arm cavities. This is so for the entire range of laser 
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FIG. 10. The 

parameters characterizing the two Fabry-Perot cavities that 
are used, together with conventional homodyne detection at 
phase 9 = n/2, to produce FD homodyne detection at the 
desired frequency-dependent phase (^). The quantities fi 
and £11 are the filters' fractional frequency offsets from the 
light's carrier frequency (|8^); Si/^y and fo/7 are the filters' 
half bandwidths in units of the half-bandwidth of the interfer- 
ometer's identical arm cavities. The functional forms of these 
parameters are Eqs. fej). 



powers, Io/Isql, that are likely to be used in QND in- 
terferometers, at least in the early years (e.g.. LIGO-III; 
ca. 2008-2010). Moreover, the filter cavities' fractional 
frequency offsets £j are of order unity (—0.5 < £j < 2). 
Thus, the desired properties of the filter cavities are not 
much different from those of the interferometer's arm cav- 
ities. 



In Sec. VI below, we shall see that the most serious 
limitation on the sensitivities of variational-output and 
squeezed-variational interferometers is optical loss in the 
filter cavities. To minimize losses, the cavities should be 
very long (so the cavities' stored light encounters the mir- 
rors a minimum number of times). This suggests placing 
the filter cavities in the interferometer's 4km-long arms, 
alongside the interferometer's arm cavities. 



C. Squeezing with frequency-dependent squeeze 
angle 

Just as the variational-output and squeezed-variational 
interferometers require homodyne detection at a FD 
phase, so a squeezed-input interferometer requires 
squeezing at a FD angle A(O). 

The nonlinear-optics techniques currently used for 
squeezing will produce a squeeze angle that is nearly 
constant over the very narrow frequency band of 
gravitational-wave interferometers, \lo — u> \ < (a few) x 
7 ~ 10 _12 w o . What we need is a way to change the 
squeeze angle from its constant nonlinear-optics-induced 
value to the desired frequency-dependent value, A = 
[Eq. ©]. 
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Just as FD homodyne detection can be achieved by 
sending the light field through appropriate filters followed 
by a frequency-independent homodyne device, so also FD 
squeezing can be achieved by squeezing the input field in 
the standard frequency-independent way, and then send- 
ing it through appropriate filters. Moreover, since the 
necessary squeeze angle has the same frequency depen- 
dence — 4>(f2) as the homodyne phase (aside from sign), 
the filters needed in FD squeezing are nearly the same 
as those needed in FD homodyne detection: The filter- 
ing can be achieved by sending the squeezed input field 
through two Fabry-Perot cavities before injecting it into 
the interferometer, and the cavity parameters are given 
by formulae analogous to Eqs. (B9h: 
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The details of the calculations are essentially the same 
as Appendix C, but with Eq. (CI) changed into the fol- 



lowing expression for the initial frequency-independent 
squeeze angle 9 and the cavities' frequency-dependent 
phase shifts aj±: 



tan$(fi) 



ft 2 ( 7 2 + Q 2 ) 
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VI. INFLUENCE OF OPTICAL LOSSES ON QND 
INTERFEROMETERS 

A. The role of losses 

It is well known that, when one is working with 
squeezed light, any source of optical loss (whether fun- 
damentally irreversible or not) can debilitate the light's 
squeezed state. This is because, wherever the squeezed 
light can leave one's optical system, vacuum field can 
(and must) enter by the inverse route; and the entering 
vacuum field will generally be unsqueezed |3l| ] . 

All of the QND interferometers discussed in this paper 
rely on squeezed-light correlations in order to beat the 



SQL — with the squeezing always produced ponderomo- 
tively inside the interferometer and, in some designs, also 
present in the dark-port input field. Thus, optical loss is 
a serious issue for all the QND interferometers. 

In this section we shall study the influence of optical 
losses on the optimized sensitivities of our three types of 
QND interferometers. 



B. Sources of optical loss 

The sources of optical loss in our interferometers are 
the following: 

• For light inside the interferometer's arm cavities 
and inside the Fabry-Perot filter cavities: scatter- 
ing and absorption on the mirrors and finite trans- 
missivity through the end mirrors. We shall discuss 
these quantitatively at the end of the present sub- 
section. (In addition, wave front errors and bire- 
fringence produced in the arm cavities and filters, 
e.g. via power-dependent changes in the shapes and 
optical properties of the mirrors, will produce mode 
miss-matching and thence losses in subsequent ele- 
ments of the output optical train.) 

• For squeezed vacuum being injected into the in- 
terferometer: fractional photon losses e c i rc in the 
circiitooT^] used to do the injection, in the beam 
splitter ebs) and in mode-matching into the inter- 
ferometer e mm . 

• For the signal light traveling out of the interferom- 
eter: In addition to losses in the arm cavities and 
filter cavities, also fractional photon losses in the 
beam splitter eb SJ in the circulator e c i rc , in mode- 
matching into each of the filter cavities e mm , in 
mode-matching with the local-oscillator light used 
in the homodyne detection ei , and in the photodi- 
ode inefficiency e p d. 

It is essential to pursue R&D with the aim of driving 
these fractional photon losses down to 



0.001 



(92) 



These loss levels are certainly daunting. However, it is 
well to keep in mind that attaining the absolute lowest 
loss levels will likely be an essential component of any 
advanced interferometer that attempts to challenge and 
surpass the SQL. In the current case, discussions with 
Stan Whitcomb and the laboratory experience of one of 



8 The circulator is a four-port optical device that separates 
spatially the injected input and the returning output from 
the interferometer; see Fig. |l| It can be implemented via a 
Faraday rotator in conjuction with two linear polarizers. 
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the authors (HJK) lead us to suggest that it may be 
technically plausible to achieve the levels of Eq. ( |9^ ) 
in the LIGO-III time frame, though a vigorous research 
effort will be needed to determine the actual feasibility 
These loss levels are certainly daunting, but based on 
discussions with Stan Whitcomb and on the laboratory 
experience of one of the authors (HJK), it seems techni- 
cally plausible that they can be achieved in the LIGO-III 
time frame. 

The arm cavities are a dangerous source of losses be- 
cause the light bounces back and forth in them so many 
times. We denote by C the probability that a photon 
in an arm cavity gets lost during one round-trip through 
the cavity, due to scattering and absorption in each of 
the two mirrors and transmission through the end mir- 
ror. With much R&D, by the LIGO-III time frame this 
loss coefficient per round trip may be as low as 



A fraction 



£ ~ 20 x 10" 



2C 



T 



2 7 L/c 



0.0012 



(93) 



(94) 



of the carrier photons that im ping e on each arm cavity 
gets lost in the cavity [cf. Eq. ( B25 ) on resonance so £ — 
e]. (Note the absence of any subscript on this particular 
e.) For si de-ba nd light the net fractional loss [denoted 
Eq. ( |100| ) below] is also of order e. 
Each filter cavity, J = I or II, has an analogous loss 
coefficient C j ~ C and fractional loss of resonant photons 



2Cj _ C 
~T7 ~ 2SjLj/c 



(95) 



Because (as we shall see), the filter cavities' losses place 
severe limits on the interferometer sensitivity, we shall 
minimize their net fractional loss in our numerical esti- 
mates by making the filter cavities as long as possible: 
L.j = L = 4km. Then the ratio of Eqs. @ and (@) 
gives 



ej = efr/6j) ~ (0.5 to 2)e 



(96) 



[cf. the last sentence of App. B; also the lossless input- 
output relation ( |l6|) and Fig. E|. Here, accurate to first 
order in e, 



/?* = arctan I - — \- ] = f3 - — £ ^ 2 — 



(98) 



is the loss-modificd| phase (3 [Eq. (|l7)], and the coupling 
coefficient is reduced slightly by the ossesf^j 



(J //sql)2 7 4 



ft 2 [ 7 2 (l + e/2) 2 + ft 2 ] 
[cf. Eq. @], where 



£ 



2 7 2 



2e 



7 2 + n 2 1 + (n/jY 



(99) 



(100) 



is the net fractio nal lo ss of sideband photons in the arm 
cavities [cf. Eq. (B25)]. Accurate to first order in the 
losses, the output quadrature noise operators in Eq. (|97| ) 
have the form 



ni 



Ab 2 = a 2 e 2 ^ (l - ±s\ + V£e^n 2 

-£*(ai + ^/7j2 ni )e 2lf3 ' (101) 

[cf. last sentence of App. B and cf. Eq. (fl6|)]. Here n± and 
n,2 are the net quadrature field amplitudes that impinge 
on the interferometer's arm cavities at their various sites 
of optical loss. We shall call rij the quadrature ampli- 
tudes of the arm cavities' loss-noise field. They are com- 
plete analogs of the input and output fields' quadrature 
amplitudes a,j and bf they are related to the loss-noise 
field's annihilation and creation operators n± and n± in 
the standard way [analog of Eqs. (g)], they have the stan- 
dard commutation relations [analog of Eqs. (j7|)], and they 
commute with the dark-port input field amplitudes aj. 

Equations (101) have a simple physical interpretation. 
The dark-port input field aj at frequency oj ± Q gets at- 
tenuated by a fractional amount £/2 while in the inter- 
ferometer (corresponding to a photon-number fractional 
loss £), and the lost field gets replaced, in the output 



C. Input-output relation for lossy interferometer 



We show in Appendix [B] that, accurate to first order in 
the arm-cavity losses (and ignoring beam-splitter losses 
which we shall deal with separately below), the relation 
between the input to the interferometer's beam splitter 
(field amplitudes aj) and the output from the beam split- 
ter (field amplitudes bj) takes the following form 

61 = A61 , 6 2 = A6 2 + y2£7-^e l/3 - (97) 

"SQL 



The loss modification, i.e. the difference between /3* and 
/3, turns out to influence the gravitational-wave noise only at 
second order in e and thus is unimportant; see Footnote |l2| 
below. 

As is discussed in Footnote [u| in Eq. @ for /C» , strictly 
speaking, I D is not the input power to the interferometer, but 
rather is the input power reduced by the losses that occur in 
the input optics, beamsplitter, and arm cavities. We ignore 
this delicacy since its only effect in our final formulas is a 
slight renormalization of 7 Q . 
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FIG. 11. The output light b is sent through a lossy out- 
put optical train, including a beam splitter, circulator, cavity 
filters I and II, a mixer with local oscillator light and a photo- 
diode. The result (aside from an unimportant phase shift a m ) 
is the desired measured quantity This actual process, 

sketched on the left side of the equality sign, is mathemati- 
cally equivalent to the idealized process sketched on the right 
side: The cavities' loss effects are introduced first, produc- 
ing b, which is then sent through an idealized, lossless optical 
train including the filters. 



light, by a small bit of loss-noise field y/£n,j. The phase 
shift [3 that the interferometer cavities put onto the loss- 
noise field is half that put onto the dark-port input field 
because of the different routes by which the aj and rij 
get into the arm cavities. 

The radiation-pressure back-action force on the test 
mass is produced by a beating of the laser's carrier light 
against the in-phase quadrature of the inside-cavity noise 
field eti + y/e/2 m . Thus, it is a\ + x/e/2 n\ that appears 
in the output light's back-action noise (last term of A6 2 ). 



D. Noise from losses in the output optical train and 
the homodyne filters 

The output quadrature operators bj get fed through an 
output optical train including the beam splitter, circula- 
tor (if present), filter cavities (if present in the output as 
opposed to the input), local-oscillator mixer, and photo- 
diode. Losses in all these elements will modify the Abj. 
In analyzing these modifications, we shall not assume, 
initially, that the FD homodyne phase is rather, 
we shall give it an arb itrary value £(fi) (as we did in our 
lossless analysis, Sec. [V C ), and shall optimize £ at the 
end. The optimal £ will turn out to be affected negligibly 
by the losses; i.e., it will still be = arccot/C. 

By analogy with the effects of arm-cavity losses [fac- 
tors £ in Eqs. (101)], the effects of the optical-train losses 
on the output fields bj can be computed in the manner 
sketched in Fig. |ll|: The process of sending the quadra- 
ture amplitudes bj through the optical train is equivalent 
to (i) sending bj through a "loss device" to obtain loss- 
modified fields bj, and then (ii) sending bj through the 
lossless optical train] 1 " 1 ] 

Because the filter cavities have frequency offsets £j 
that make their losses different in the upper and lower 



n Yanbei Chen has shown that it does not matter 

whether the losses are placed before or after the lossless train. 



side bands, the influence of the losses is most simply ex- 
pressed in terms of the annihilation operators for the side 
bands b± , rather than in terms of the quadrature ampli- 
tudes bj. In terms of b±, the equation describing the 
influence of losses is identical to that in the case of the 



arm cavities with fixed mirrors, Eqs. (101) with JC = 



(102) 



Here (i) the sum is over the two filter cavities J = I and II 
(which must be treated specially) and over the rest of the 
output optical train, denoted J = OT; (ii) £j is the net 
fractional loss of photons in element J; (iii) nj± is the 
annihilation operator for the loss-noise field introduced 
by element J; (iv) for each filter I or II, the analog of the 
phase factor 2(3 of Eq. ( |101| ) gets put onto the light in 
the subsequent lossless filter and thus is absent here; and 
(v) we have absorbed a phase factor into the definition 
of n J± . 

The net fractional photon loss in a fil ter c avity must 
be identical to that in an arm cavity, Eq. (100), if written 
in terms of the cavity's half bandwidth (7 for arm cavity, 
5 j for filter cavity) and the difference between the field's 
frequency to = uj ±Q and the cavity's resonant frequency 



u> Tes {bj — U! Ies = ±f2 for arm cavity 



for filter cavity). Therefore, Eq. (10C) implies that 



ijSj ± n 



2e. ; 



-j± 



1 + (o ± n/5j 



for J = I, II 



(103) 



For the remainder of the optical train, the net fractional 
photon loss Sot is the sum of the contributions from the 
various elements and is independent of frequency: 

£0T± = £oT = £bs + Ecirc + 2e mm + Cfo + £ p d ~ 0.006 . 

(104) 

By expressing b± and nj± in terms of bj and njj (for 
j = 1,2) via the analog of Eq. (||), inserting these expres- 
sions into Eq. (102), then computing bj via the analog of 
Eq. (pi), we obtain 



h = \l- 5 Wj h-- J2(£j+ ~ Sj-)bi 



(105a) 



n J2 



(105b) 



1 J 



Here 
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^OTF = ~ y^(^J+ + 



Sot + 2(^1+ + ft- + &+ + Si-) 



^ot + ^E E i 



J=I,II s=J 



(106) 



is the net, ^-dependent loss factor for the entire output 
optic al trai n in cluding the filter cavities. From Eqs. (94), 
(EM and <KM) and Fig. pi, we infer that 



TF 



0.009 



(107) 



with only a weak dependence on frequency, which we 
shall neglect. 



In Eqs. ( 105 ), the terms i x (quantity linear in £j±)bj 
[the 62 term in 61 and the &i term in 6 2 ] wm contribute 
amounts second order in the losses (oc £j) to the signal 
and/or noise, and thus can be neglected. We shall flag 
our neglect of these terms below, when they arise. 



E. Computation of noise spectra for 
variational-output and squeezed-variational 
interferometers 

The output of a squeezed-variational interferometer 
or variational-output interferometer is the frequency- 
dependent quadrature b( depicted in Fig. [n]. This quan- 
tity, when split into signal oc h plus noise oc A6^, takes 
the following form: 



b^ = 61 cos £ + 62 sin ( 



sin£ 



\ F \ h 

2 / /J S Q L 



,03, 



+ 



sin C 



(108) 



cf. Eqs. @ and ( [ToB] ). Here we have omitted an imag- 
inary part of the factor 1 — ^£qtf [arising from the 62 




term in 61, Eq. (105a)] because its modulus is second 
order in the losses (oc £j) and therefore it contributes 
negligibly to the s ignal strength. 

Equation (108) implies that the gravitational- wave 
noise operator is 

(Ab 2 +A6icotC) , 
(109) 

where we have used Eq. (^) for /C*. 

For a squeezed-variational interferometer, the dark- 
port input field a.j is in a squeezed state, with squeeze 
factor R and squeeze angle A(tt) (which, after optimiza- 
tion, will turn out to be A = ir/2 as for a lossless inter- 
ferometer). For a variational-output interferometer, aj is 
in its vacuum state, which corresponds to squeezing with 
R = so we lose no generality by assuming a squeezed 



input. Since all the noise fields except a,j are in their 
vacuum states, the light's full input state is 



= |0„) ® |0„ OT ) <g> |0 ni ) <g> |0 niI ) (» S(R, A)|0 8 ) 



where the notation should be obvious. 

The gravitational-wave noise is proportional to 

(m|/i n /i„'|m) = (0\h ns h ns '\0) 



(110) 



(111) 



where |0) is the vacuum state of all the noise fields a, 
n, noT) n\, and nn; and h ns is the usual squeezed noise 
operator 



S*(R, X)h n S(R, A) 



(112) 



We bring this squeezed-n oise operator i nto an explicit 
form by (i) inserting Eq. (109) into Eq. ( |112| ), then (ii) 
replacing the A&j's by expressions ( |l05|) [with A put onto 
all the 6's, i.e. with the signal removed] , then (iii) replac- 
ing th e A frj's by expressions ( |l0lD , and then (iv) invoking 
Eqs. (A8) for the action of the squeeze operators on the 



aj's. The result is 



K. s = ( 1 - -£ ) x [Eq. (69)] 



{ (-/Ce 4/3 + VS cot c) e i0 m + Ve, 

E + v/^)cotC- i ( v / ^7- \f£~T-) 



nji 




where we have omitted terms , aris ing from 62 in Eq. 
(|l05a|) and from h in Eq. ( |l05b|) , which contribute 



amounts 0(£j) to S^; and we have omitted a termP] pro- 
portional to /?* — (3 which contributes an amount C(e 2 ). 

By virtue of Eq. (Ill) and the argument preceding Eqs. 
( ptf ) , we can regard all of the quadrature noise operators 
a,j, rij, njj in this h ns as random processes with unit 
spectral densities and vanishing cross-spectral densities. 
Correspondingly, the gravitational-wave noise is the sum 
of the squared moduli o f the coefficients of the quadrature 
noise operators in Eq. (113): 



12 this term is an imaginary part, 2i(/3„ — f3)K, = 
— |ie/Csin2/3, of the quantity K,, which enters Eq. (p9) via 
Eq. ([70)). Because this imaginary part produces a correction 
to the loss-free part of h n that is 90 degrees out of phase with 
the loss-free part and is of order e, it produces a correction to 
Sh that is quadratic in e and thus negligible. 
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2JC 

x{e- 2R 



sinh2i?[l-cos2($ + A)]} (114) 



+/C 2 - + (1 - K cot QE + (1 + cot 2 QEotf 



where 



JC = JC — cot C , $ = arccotA^ 



(115) 



[Eq. (|70|)1 . In Eq. (114), the first two lines come from a\ 
and ai [squeezed vacuum entering the dark port; cf. Eq. 
(fn])] modified by losses in the arm cavities [the factor 
1 — £/2)]; the first two terms on the third line come from 
ni and n-i [shot noise due to vacuum entering at loss 
points in the arm cavities] ; and the last term comes from 
nji and nj2 [shot noise due to vacuum entering at loss 
points in the output optical train, including the filters] . 

As for the lossless interferometer [Eqs. ( [72] ) and (f73f)], 
the noise ( 114 ) is minimized by setting the input squeeze 
angle A and output homodyne phase £ to 



A = 7r/2 , ( = $ = arccot/C 



(116) 



[aside from a neglible correction <5£ = (£ + 2£oTp)e~ 2R / 
QC + /C -1 )]. This optimization produces K — and A = 

$ = 7T/2, so 



S h = 



h 2 

"SQL 



(117) 



(i - ¥) 



-2R 



Note that the optimization has entailed a squeezed in- 
put with frequency-independent squeeze phase, as in the 
lossless interferometer; so no filters are needed in the in- 
put. The output filters must produce a FD homodyne 
angle ( = that is the same as in the lossless case 

and therefore can be achieved by two long, Fabry-Perot 
cavities. 



It is instructive to compare the noise ( 117 ) for a lossy 
squeezed- variational interferometer with that ( f73| ) for one 
without optical losses. In the absence of losses, the out- 
put's FD homodyne detection can completely remove the 
radiation-pressure back-action noise from the signal; only 
the shot noise, cx 1/JC oc l// , remains. Losses in the 
interferometer's arm mirrors prevent this back-action re- 
moval from being perfect: they enable a bit of vacuum 
field n to leak into the arm cavities, and this field pro- 
duces radiation-pressure noise that remains in the output 
after the FD homodyne detection (the ICe/2 term in Eq. 



The ZC£otf noise in Eq. (117) has the same dependence 
on laser power, cx JC oc I , as the radiation-pressure noise. 
Nevertheless, it is actually shot noise, not radiation pres- 
sure noise. It is produced by the vacuum loss-noise fields 
that leak into the output signal light when it encounters 
each lossy optical element. Those fields' shot noise gets 



weighted by the factor cot C = cot $ = K in the homo- 
dyne process, which accounts for their proportionality to 
fC cx I Q . 

A reasonable estimate for the amount of input-light 
squeezing that might be achieved in LIGO-III is |3j| 

e~ 2R ~ 0.1 . (118) 

By contrast, Eqs. @, @, p0| ) and ( |l06| ) suggest 

(8 + £ tf) ~ 0.01 . (119) 

This motivates our neglecting £+£otf comp ared to e~ 2R 
in expression (jH7|), and rewriting the noise (117) as 



h 2 

""SQL 



,-2R 



K. 



where 



£otf - 0.0010 



(120) 



(121) 



cf. Eqs. @ and p7| ). 

Equation ( |l20| ) is our final form for the noise spectrum 
of a lossy squeezed-variational interferometer. When we 
set the input squeeze factor to unity e~ 2R = 1, it be- 
comes the noise spectrum for a lossy variational-output 
interferometer: 



S h 



h 2 

"SQL 



1 

K 



(122) 



Errors A A = A — tt/2 in the input squeeze angle and 
AC = £ — arccot/C in the output homodyne phase will in- 
crease the noise spectral density. By performing a power 
series expansion of expression (114), we obtain for the 
noise increase 

h 2 

"SQL 



AS h = [sinh 2R AX 2 - 2(1 + K?) sinh 2R AAA( 

(l + lC) 2 e 2R 



AC 



(123) 



^e 2 «[AA-(l+/C 2 )AC] 2 , 



where the second expression is accurate in the limit 



-2R 



e'^ 3> e~'"\ Numerical evaluations show that, for 
e~ 2R — 0.1 and e* = 0.01 (see above), an d fo r K ~ 1 
to 3 (the range of greatest interest; cf. Sec. VII), Ay/Sh 
will be less than \\f~S~h so long as: (i) the input squeeze 
angle is accurate to |AA| < 0.05, and (ii) the FD out- 
put homodyne phase is accurate to \AC\ < 0.01. At 
JC = 1 the FD phase's required accuracy is reduced to 
|AC| < 0.04. The FD phase C is determined by the filter 
cavities' half bandwidths Sj and fractional frequency off- 
sets £j, and the local oscillator phase or equivalently the 
final, conventional homodyne detector's homodyne phase 
9. The filter cavities' half bandwidths Sj (or equivalently 
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their finesses) are fixed by the mirror coatings. Coating- 
produced errors in 8j can be compensated to some de- 
gree by tuning the fractional frequency offsets £j (via 
adjusting the mirror positions) and by tuning the local 
oscillator phase or equivalently 9. Finesse errors as large 
as five per cent, \A5j\/Sj < 0.05, can be compensated to 
yield the required \A(\ < 0.01 by tuning the offsets and 
homodyne phase to one perce nt a ccuracy, |A£j| < 0.01, 
A0 < 0.01 [Eqs. (g), @, (|C2|) and Fig. @. These 
requirements are challenging. 



F. Computation of the noise spectrum for a 
squeezed-input interferometer 

For a squeezed-input interferometer, as for squeezed- 
variational, the losses in the input optical train influence 
the noise only through their impact on the squeeze factor 
e~ 2R nu 0.1 of the dark-port vacuum when it enters the 
arm cavities. By contrast, losses in the arm cavities and 
in the output optical train will produce noise in much 
the same manner as they do for a squeezed-variational 
interferometer. More specifically: 

The effect of arm-cavity and output-train losses on 
the squeezed noise operator h ns c an be read off of the 
squeezed-variational formula (113) as follows: (i) Set 
£ = 7r/2 so the quantity measured is 62 [no output filter- 
ing; Eq. ( |108] )]; (ii) correspondingly set cot £ = 0, K = /C, 
and $ = $ = arccot/C [Eqs. ( |ll5| )1; (iii) in the sum over 



J include only J = OT and not J = I, II since there are 
no output filters. The result is 



[Eq. (44)] 




n2 



V £oT' fiOT'2 

(124) 



Here the prime on the subscript OT indicates that we 
must omit losses due to mode-matching into the output 
filters and mixing with the local oscillator, since there are 
no output filters or homodyne detection. Correspond- 
ingly, 



£ OT' = £bs 



0.003 



(125) 



is the net fractional photon loss in the output optical 
train. 

Treating the quadrature noise operators as random 
processes with unit spectral density and va nishi ng cross 
spectral densities, we read off Sh from Eq. (124): 



Sh 



h 2 

"■SQL 



£ + £0T' 

K 



x {cosh2i? - cos[2(A + $)] sinh 2R} 



(126) 



As in the lossles case, the noise is minimized by squeez- 
ing the dark-port input at the FD angle A(f2) = — $ = 
-arccot/C [Eq. @]. The result is 
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(127) 



For our estimated squeezing e ~ 0.1 and losses 
£ OT , ~ £ ~ e < 0.003 in the LIGO-III time frame, the 
loss parameters are small compared to the squeeze, and 
thus contribute negligibly to the noise, so Sh is well ap- 
proximated by the lossless formula 
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(128) 



However, it is important to keep in mind that the input 
squeeze factor e~ 2R is constrained not only by the physics 
of the squeezing appratus, but also by losses in the input 
optical train and mode matc hing into the arm cavities. 

By expanding expression ( 126 ) in powers of A A = A + 
arccot/C, we see that the fractional increase in noise due 
to errors in the FD squeeze angle is 



^Vs~ h 
VSh~ 



e 2R sinh 2R SX 2 



./Hi 



-AX 



(129) 



For e =0.1, this fractional noise increase will be less 
than 1/4 so long as AA is less than 0.07. This translates 
into accuracies of ~ 7 per cent for the prefilter squeeze 
angle, ~ 15 per cent for the filter cavities' fractional fre- 
quency offsets (|A£j| < 0.15), and ~ 10 per cent for 
the cavities' half bandwidths or equivalently their finesses 
(ASj/Sj < 0.1). These constraints are significantly less 
severe than thos e for a squeezed-variational interferom- 
eter (end of Sec. VIE ); but, as we shall see, the poten- 
tial performance of this squeezed-input interferometer is 
poorer by a factor ~ 1.5-2 than that of the squeezed- 
variational one. 



VII. DISCUSSION OF THE 
INTERFEROMETERS' NOISE SPECTRA 

The noise spectra for our three lossy QND interferom- 



eters, Eqs. (120), (122) and (128), all have the same uni- 
versal form — a form identical to that for a conventional 
broadband interferometer, Eq. (27). Only the parameters 
fi and a characterizing the noise differ from one interfer- 
ometer to another. This universal form can be written 
as 



+ 



(130a) 



where f2* is the following function of angular frequency 
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FIG. 12. Universal nois e curve for conventional and QND 
interferometers [Eqs. ( 130 )]. 
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FIG. 13. Universal noise curve plotted as a function of an- 
gular frequency Q for various values of the dimensionless fre- 
quency parameter a. 
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(130b) 



and 1isql(Q) is given by Eq. (|20j). Notice that f2» = 1 
when Q = 7 ~ 100 Hz; fi, = (n/-y)/y/2 when fi < 7, 
and Q* = (fi/7) 2 /v / 2 when > 7. 

This universal noise curve is plotted as a function of 
£1* in Fig. |l2| Its two parameters are the minimum value 
/i of the noise, i.e. the minimum amplitude noise in units 
of the SQL, and the dimensionless frequency a (in units 
of f2») at which the noise takes on this minimum value. 

Figure [l3| shows this universal noise curve plotted as 
a function of an gular frequency O. Notice that, because 
of the relation ( 130b ) between fi* and f2, the shape of 
the noise curve depends modestly on the location a of its 
minimum. 

The values of the parameters [i and a for our various 
interferometer configurations are shown in Table y. No- 
tice the following details of this table: (i) The minimum 
noise /i (the optimal amount by which the SQL can be 
beat) is independent of the laser input power I Q in all 
cases; it depends only on the level of input squeezing 
e~ 2R and the level of losses e*. (ii) For our estimated 
loss level and squeeze level, the squeezed-input inter- 
ferometer and variational-output interferometer achieve 
the same [i ~ 0.32, while the squeezed-variational inter- 



TABLE II. The values of the parameters /1 = (minimum 
noise) and a — (frequency of minimum) for various interfer- 
ometer ("IFO") configurations: Conv = C onve ntional broad- 
band [Eq. @], SI = Squeezed-Input [Eq. (fl2S|)], VO = Vari- 
ation al- Ou tput [Eq. (122)], and SV = Squeezed- Variational 
[Eq. (120)]. The numerical values are for e~ 2R = 0.1 and 
e* = 0.01. 



IFO 

Conv. 
SI 

VO 
SV 



1 



Ve~ 2R ~ 0.32 



y Io/Isql 



1/4 



0.32 



(e- 2fl e,) 1/4 



0.18 



W-fSQL 
W-fSQL 



ye. 



i/ f SQL 



3.2 



ferometer achieves a moderately lower /j, ~ 0.18. (iii) 
The frequency f2» = a at which the minimum noise 
is achieved is proportional to */ I /Isql- (Recall that 
7gQL is the input power required for a conventional in- 
terferometer to reach the SQL at the angular frequency 
n = 7 ~ 2ir x 100 Hz, i.e. at Q* = 1; to do so, the 
conventional interferometer must have (7=1.) (iv) For 
lo = /sql, the squeezed-input interferometer has a = 1, 
but the variational-output and squeezed-variational in- 
terferometers have a < 1, which means that the mini- 
mum of the noise curve is at Q < 7 ~ 100 Hz. To push a 
up to unity, i.e. to push the noise-curve minimum up to 
Q = 7, requires I /Isql — — 10 in a variational- 

output interferometer, and I /Isql = \J e~ 2R ~Je* — 3.2 
in a squeezed-variational interferometer. 

The importance of pushing a up to unity or higher is 
explained in Fig. 
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This figure requires some discussion: 
The most promising gravitational waves for LIGO are 
those from the last few minutes of inspiral of black- 
hole/ black- hole binaries, black- hole/neutron-star bina- 
ries, and neutron-star/neutron-star binaries. The ampli- 
tude signal-to- noise ratio S/N produced by these waves 
is given by 



N 2 



= 4 



\h\ 2 dfl _ ] 



\nh\ 2 dlntt 



(131) 



where h is the Fourier transform of the waveform h(t). 
For the inspiraling binary \Qh\ is nearly independent of 
frequency throughout the LIGO band 0] , so the signal- 
to-noise ratio is optimized by making QSh(&) as small as 
possible over as wide a range of In as possible. 

Figure H plots y/ns h (fi.) as a function of fi/7 using 
logarithmic scales on both axes, and using the minimum- 
noise parameter fi = 0.18 corresponding to our fiducial 
squeezed-variational interferometer (though the specific 
value of /i is irrelevant to our present discussion). From 
the shapes of the curves it should be evident that the 
larger is the frequency of the noise minimum, i.e. the 
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FIG. 14. Noise curves for SQL interferometers with noise 
minima fi = 0.18 and various values of the freq uency param- 
eter a. The vertical axis is weighted by so the curves 
give an indication of the relative noise in searches for waves 
from inspiraling binaries; see text. The noise curves are la- 
beled by the power I /Isql required by a squeezed- variational 
interferometer to achieve the given a. 



larger is a at fixed fi, the larger will be the S/N for in- 
spiraling binaries. 

A second factor dictates using large <j, in particular 
er > 1. This is thermal noise in the interferometer's test- 
mass suspension fibers. The thermal noise scales with fre- 



quency as y / fiS'* hermal (tt) oc fT 2 or oc fl~ 5 / 2 depending 

on the nature of the dissipation J35| ; see the steep dashed 
curve in Fig. |lj. It seems realistic to expect, in LIGO- 
III, that this thermal noise will be at approximately the 
level shown in the figure, so it compromises the perfor- 
mance of QND interferometers at f2 < 0.57 ~ 50 Hz 
|t|,[36|. Correspondingly, to avoid the thermal noise sig- 
nificantly debilitating the S/N for inspiraling binaries, it 
will be necessary to have cr > 1. 

Because a scales as \J I /Isql for all interferometer 
designs, large a entails large laser power. In particular, 
a > 1 requires I > /sql! cf. Table Q. For our fiducial 
parameters (Table Q), Isql = 10 kW, which corresponds 
to an optical power circulating in each of the interferom- 
eter's arm cavities 



W. 



SQL _ -fsQL/2 _ mcLj 3 



8u> 



= 0.62 MW 



(132) 



To construct mirrors capable of handling this huge power 
will be an enormous technical challenge (even though 
this is approximately the circulating power contemplated 
for LIGO-II). To operate with a circulating power much 
larger than this might not be possible. Therefore, it may 
be important in LIGO-III to achieve a > 1 while keeping 
Io/Isql not much larger than unity. 

The squeezed-input interferometer, with its a = 
y/Io/IsQL (Table |0]) is the most attractive from 
this point of view [and also in terms of its re- 
quire d filt er and squeeze-phase accuracies; cf. end of 
Sec. VTF]; and the variational-output with its a = 

y/ ^Io/Isql — \/0.1/ o /Jsql is the least attrac- 
tive. The squeezed-variational interferometer, with a = 



\fy/e*/e ^Io/Isql — \/0-S2I o /I S q L requires a mod- 
estly higher laser power to reach a = 1 than the squeezed- 
input [and requires better filter and squeeze-phase ac- 
curacies], but it is capable of a lower noise minimum, 
fx ~ (e-^e*) 1 / 4 ~ 0.18 vs. fi = ^e~ 2R ~ 0.32 for 
squeezed-input. 

This suggests an R&D strategy: Focus on input 
squeezing as a key foundation for LIGO-III (it is needed 
both for squeezed-input and squeezed-variational inter- 
ferometers), and in parallel (i) develop the technology 
and techniques for the FD homodyne detection required 
by squeezed-variational configurations, (ii) work to drive 
down optical losses to the levels e ~ e c i rc ~ £bs ~ £mm ~ 
eio ~ e P d ~ 0.001 [Eq. (p2])j, and (since ponderomotive 
squeezing, which underlies all our QND interferometers, 
has never been seen) (iii) carry out experiments in a small 
test appratus to demonstrate ponderomotive squeezing 
and to search for unexpected obstacles and imperfections 
in it. 

If both input squeezing and FD homodyne detec- 
tion can be implemented successfully, then the squeezed- 
variational interferometer is likely to achieve better per- 
formance than any other configuration discussed in this 
paper, despite its apparent need for higher laser power 
(e.g. Io/Isql — 3.2 to achieve a = 1 compared to 
Io/Isql = 1 for squeezed input, with our fiducial pa- 
rameters). If powers as high as I /Isql — 3.2 cannot be 
handled, then we can operate the squeezed-variational 
interferometer with a lower power without much loss of 
performance. 

Consider, for example, y/Sh/hsQh evaluated at fl — 
o~ — 2,tt x 100 Hz, as a function of I /Isql in a variational- 
output interferometer with our fiducial e~ 2R = 0.1 
and e„ = 0.01. The optimal I /Isql — 3.2 produces 
\/Sh{j)/h S Q L = 0.18; pushing 7 //sql down by a factor 
2, to 1.6, increases the noise at f2 = 7 by only 10 per- 
cent, to 0.20; pushing down all the way to I /IsQh 
1 increases the noise to only V^XtJ/^sql = 0.23, 
which is still significantly lower noise than the optimized 
squeezed-input interferometer (0.32 at I /Isql = 1)- 

It is worth recalling that for noncosmological sources 
(sources at distance 3 Gpc), the volume of the universe 
that can be searched for a given type of source scales as 
the inverse cube of the amplitude noise, so a noise level 
^/Sh/h-SQL — 0.18 corresponds to search- volume increase 
of 1/0. 18 3 ~ 180 over a SQL- limited interferometer, i.e. 
over LIGO-II. 



VIII. CONCLUSIONS 

In this paper we have explored three candidate 
ideas for QND LIGO-III interferometers: squeezed- 
input, variational-output, and squeezed variational. The 
squeezed-input and squeezed-variational interferometers 
both looks quite promising. For our estimated levels of 
optical loss and levels of squeezing, and for an input laser 
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power I /Isql = 1 (the LIGO-II level), the squeezed- 
input interferometer could achieve a noise fi ~ 0.32 of the 
SQL, with a corresponding increase V ~ 1/0. 32 3 ~ 30 
over LIGO-II in the volume of the universe that could 
be searched for a given source, at non-cosmological dis- 
tances. The squeezed-variational interferometer could 
achieve /i ~ 0.23 of the SQL with a search-volume in- 
crease over LIGO-II of V ~ 80. If the optics can handle a 
laser power I /IsQh — 3.2, then the squeezed-variational 
interferometer could reach fi ~ 0.18 of the SQL and a 
search- volume increase of V ~ 180. These numbers scale 
with the losses, squeezing, and laser power as shown in 
Table [fi[ 

The squeezed-input and squeezed-variational designs 
are therefore sufficiently promising to merit serious fur- 
ther study. Some of the issues that need theoretical anal- 
ysis are: 

• How can one incorporate into these interferometer 
designs the various light modulations that are re- 
quired, in a real gravitational-wave interferometer, 
to (i) make the interferometer be shot-noise lim- 
ited (put the gravitational- wave signal into ^100 
Hz sidebands of a MHz modulation]^), (ii) control 
the mirror positions and orientations, ... |37|| . 

• What accuracies and other characteristics are 
needed for the interferometers' new elements: the 
circulator, filter cavities,]^] and input squeezing? 
How can these be achieved? For example, how sta- 
ble must be the local oscillator for the conventional 
homodyne detector, and can it be achieved simply 
by tapping some light off the interferometer's out- 
put or input beam? 

• If the filter cavities are placed in the same long 
vacuum tubes as the interferometer's arm cavities 
(with their enormous circulating power), what will 
be the nature and level of noise due to scattering 
of light from the test-mass cavites to the filter cav- 
ities? (We thank Eanna Flanagan for raising this 
issue.) 

• Can the filter cavities be made to serve multiple 
purposes? For example, is it possible to use a single 



optical cavity for both filters, e.g., with the two 
filters based on two different polarization states (for 
which the filter might be made to behave differently 
via birefringence), or with the two filters based on 
different, adjacent longitudinal modes? As another 
example, could an output filter cavity be used as 
a source of pondcromotively squeezed vacuum for 
input into the interferometer's dark port?P| 

• Signal recycling via resonant-sideband extraction 
(RSE) H is likely to be a standard tool in LIGO- 
II Q. How can one best implement RSE simulta- 
neously with the FD homodyne detection (and in- 
put squeezing) of a variational-output (or squeezed- 
variational) interferometer? |37j How can one best 
achieve the FD homodyne's filtration [which will 
entail a different frequency dependence from 
that in this paper's non-RSE designs]? 



In this paper's analysis we have made a number of 
simplifying approximations [e.g., our approximat- 
ing the phase of the coefficient of fj in Eq. (B24) 
by 2/3 an approximation that fails by a frequency- 
dependent amount which can be nearly as large as 
one per cent] . At what level of sensitivity do these 
approximations become problematic (e.g., for our 
proposed two-cavity way of achieving the necessary 
FD homodyne detection), and how can the result- 
ing problems be overcome? 



• Our analysis is based on the crucial assumption 
that the interferometer's output is strictly linear 
in its input g. Matsko and Vyatchanin [g^| have 
shown that this is not quite correct. In the in- 
terferometer's arms the back-action-induced mir- 
ror displacement X produces a phase shift of re- 
flected light given by e - 2tnx / c ^ which our lin- 
earized analysis approximates as 1 — 2iQX/c [cf. 
Eq. ( B10| )]; when the better approximation 1 — 
2iflX/c— 2(fLY/c) 2 is used, the result is additional, 
nonlinear noise, which limits the cancellation of the 
back-action noise by the shot noise and produces a 
limit M 



h 2 

"■SQL 



2 jVc 



1/5 



5 x 10 



-5/2 
"SQL 



(133) 



SQL 



13 LIGO scientists are currently exploring the possibil- 
ity of achieving shot-noise-limited performance in LIGO- 
II without this modulation/demodulation. The modula- 
tion/demodulation may, in fact, be replaced in LIGO-II by 
homodyne detection at the interferometer output, making it 
more nearly like our paper's LIGO-III designs. 

14 The filter cavities will require a mechanical stability far 
less than that of the arm cavities, since the carrier power in 
the output light is small and filter mirror displacements of 
magnitude ~ hL therefore do not imprint a significant signal 
on the light. 



on the sensitivity that any of our QND designs can 
achieve. Here 



15 For ponderomotively squeezed vacuum, the squeeze angle 
is frequency dependent, with d(j)/dQ of the opposite sign to 
that needed by a squeezed-input interferometer. This must 
be compensated by a filtering different from that discussed in 
Sec. 0. 
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is the number of quanta entering a SQL interfer- 
omet er in time 7 -1 ~ 2 ms. The nonlinear limita- 
tion (133) is sufficiently far below the SQL that 
we need not be concerned about it. Are there 
any other, more serious sources of nonlinearity that 
might compromise the performance of these inter- 
ferometers? 

Experimental studies are also needed as foundations 
for any possible implementation of variational-output or 
squeezed- variational interferometers |4C[ |. Examples are 

• Studies of the debilitating effects of very high circu- 
lating powers, W ~ a few MW, and how to control 
them. 

• A continuation of efforts to achieve large squeez- 
ing, robustly, via nonlinear optics [j33| , and explo- 
ration of-|thc possibility to do so ponderomotively 

• A continuation of efforts to achieve low levels of 
losses in optical cavities and interferometers, so as 
to minimize the contamination of squeezed light by 
ordinary vacuum p5| . 

• Prototyping of FD homodyne detection by the tech- 
nique proposed in this paper: filtration followed by 
conventional homodyne detection. 

In the meantime, and in parallel with such studies, it 
is important to push hard on the effort to find practical 
QND designs that entail circulating light powers well be- 
low 1 MW jl3| , and that might be much less constrained 
by optical losses than the designs explored in this paper. 
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APPENDIX A: ROTATION AND SQUEEZE 
OPERATORS 

In this paper we make extensive use of squeeze opera- 
tors and some use of rotation operators. In this appendix 
we list properties of these operators that are be useful in 
verifying statements made the text. This appendix is 
based on the formalism for 2-photon quantum optics de- 
veloped by Caves and Schumaker p^ , |27| . 

The rotation operator R(0), which acts on the Hilbcrt 
space of the modes with frequencies uj = u> ±fl, is defined 
by 



R(9) = exp[-i6(a\a + + a [ _a-)] 



(Al) 



(Eq. (4.33) of f2q| ); here a± are the annihilation opera- 
tors, and a± the creation operators for photons in these 
modes. This operator is unitary and has the inverse 

R~ 1 (9) = R^(6) = R(—0) . (A2) 

The effect of a rotation on the modes' annihilation oper- 
ators is 



R(9)a ± R' i (9) = a ± e ie 



(A3) 



(Eq. (4.35) of p6[), and its effect on the two-photon 
quadrature amplitudes [Eqs. (0)] is 



R{Q)a 1 R\6) — a\ cos 9 — a,2 sin 9 
R(9)a 2 R i ' (9) = 01 sin 9 + a 2 cos 9 



(A4) 



(Eq. (4.36) of f§). 

The squeeze operator also acts on the Hilbert space of 
modes with frequencies — c<j ± O, and is defined by 

S(r, 4) = exp[r(a + a_e- 2 ^ - a^ate**)] (A5) 

(Eq. (4.9) of H|; Eq. (1.8) of Q). This squeeze operator 
is unitary and its inverse is 

fif-^r, <j>) = S+(r, (f>) = S(-r, <f>) = S(r, 4> + tt/2) (A6) 



(Eq. (1.9) of |27|]). The effect of a squeeze on the modes' 
annihilation operators is 

S(r,(f>)a±S^r,(t>) = a± coshr + a^e 2 ^ sinhr (A7) 

(Eq. (4.10) of pH|). From this equation and the definition 
(^) of the quadrature amplitudes, we infer the effect of a 
squeeze on those amplitudes 

S(r, <f))aiS'(r, 4>) — a\ (coshr + sinh r cos 2(f)) 
+C12 sinhr sin 2(f) , 

S(r, 0)a2>S^(r, ip) = a 2 (coshr — sinh r cos 2(j>) 

+ai sinh r sin 2(f) . (A8) 
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APPENDIX B: INPUT-OUTPUT RELATIONS 
FOR INTERFEROMETERS 

In this appendix we shall derive the input-output re- 
lations for the fields a,j and bj that enter and leave the 
interferometer's dark port. From the outset we shall in- 
clude optical losses in our derivation, the reby obtaining 
the lossy input-output relations (|97| ) and (101); the loss- 
less input-output relations (|l6|) then follow by setting 
e = 0. 



1. Fields at beam splitter 

We describe the field amplitudes entering and leaving 
the beam splitter by the notation shown in Fig. [l5| (cf. 
Fig. ||. We idealize the beam splitter as lossless in this 
appendix, and deal with its losses in the body of the 
paper in the manner sketched in Fig. [ll]. The ampli- 
tudes D & d of the field entering the beam splitter from 
the laser arc defined by the following formulae for the 
positive-frequency part of the electric field 



E; 
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[cf. Eq. 



E„ 
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and for the total electric field 
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(B2) 



Thus, D is the classical amplitude of the laser light (car- 
rier with frequency u> ), d± are the annihilation operators 
for the uj ± f2 sidebands, and d\ and c?2 are the quadra- 
ture amplitudes for the side bands. [Notice that the fac- 
tor out front is a \/2tt in Eq. (Bl) but in (p32j), and 
notice the \[2D in (B2).] The light power I a impinging 



on the beam splitter is related to the classical amplitude 
Dby 



4tt ' 



■Ac = hu D 2 



(B3) 



where the overbar means time average. (Note that D 2 
has dimensions Hz = 1/sec.) 

For all other fields the classical amplitude and sideband 
amplitudes are as indicated in the figure; for example, the 
field going toward the east cavity has classical amplitude 
D/^/2 and quadrature amplitudes /f, /f . 

With an appropriate choice of conventions fig ], the 
fields' junction conditions at the splitter are 
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FIG. 15. Field amplitudes entering and leaving the beam 
splitter (which here is idealized as lossless). The various am- 
plitudes are denned in Eqs. (B1)-(B4). 
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(B4) 



Here j — 1 or 2. 



2. Arm cavities and fields 

The east and north arm cavities are presumed to be 
identical, with power reflection and transmission coeffi- 
cients R and T for the front mirror, and R and T for the 
back mirror. The amplitude reflection and transmission 
coefficients are chosen be real, with signs — 

{+VT,-\/fi} for light that impinges on a mirror from 

outside the cavity; and {+VT, +^R}, {+ VT, +\/r} for 
light that impinges from inside the cavity. 

The dominant optical losses are for l ight im pinging on 
mirrors from inside the cavity (cf. Sec. VI B ). The influ- 
ence of the losses on the interferometer's signal and noise 
are independent of the physical nature of the losses — 
whether it is light scattering off a mirror, absorption in 
the mirror, or transmission through the end mirror. (We 
ignore the effects of mirror heating.) For computational 
simplicity, we model all the losses as due to finite trans- 
missivity C = T ^ of the end mirror, and correspond- 
ingly we set 



R + T =1 , R + T = 1 



(B5) 



The fractional loss of photons in each round trip in the 
cavity is then T, and the net fractional loss of photons in 
the arm cavities is 



2C 



2T 



(B6) 



cf. Eqs. (||) and (g). Recall that T ~ 0.033 and e 
0.0012, and also that Q ~ 7 = Tc/AL [Eqs. (||), (|| 
correspondingly, we shall make the approximations 
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FIG. 16. Field amplitudes entering and leaving an arm cav- 
ity. The cavity's front-port input and outp ut amplitudes D, 
f and g are denned in Eqs. (B3) and (B4) and Fig. |l5|, and 
it's back-port input q is defined in Eq. (Bll). 



f < T = 4 7 L/c - flL/c < 1 



(B7) 



throughout our analysis. 

Figure [l6] shows an arm cavity and the amplitudes of 
the fields that impinge on or depart from its mirrors. The 
amplitudes are those at the (front or back) mirror loca- 
tion, and the mirrors, like the beam splitter, are idealized 
as infinitesimalfy thin. 

For pedagogical simplicity, the distance from the beam 
splitter to the front mirror of each arm cavity is set to 
an integral multiple of the carrier wavelength and is as- 
sumed to be far smaller than c/Q (the wavelength asso- 
ciated with the sidebands). This means that there are 
no net phase shifts of the light in traveling between the 
beam splitter and the cavity's front mirror; i.e., the field 
amplitudes D / ' \[2 h / (or Dj\[2 & g) arriving at (or de- 
parting from) the mirror are the same as those departing 
from (or arriving at) the beam splitter; cf. Figs. 1(] and 

n 

The cavity's length is adjusted to an integral num- 
ber of carrier wavelengths so there is no carrier phase 
shift from one end of the cavity to the other, and in- 
side the cavity the carrier amplitude is amplified by 
the standard resonance factor 2/VT. (L osses are small 
enough to be of little importance for the carrier.) Be- 
cause the side bands inside the cavity have a frequency 
dependence gje~ lflt cos(u) t) at the front mirror loca- 
tion [cf. Eq. (]B2j)], they propagate down the cavity as 
Sj e ~ 4f ^ t_z ' cos{u) (t — z)] and upon reaching the back mir- 
ror (where cos[uj (t — z)] = cos[u>i\), they have acquired 
the phase shift indicated in the figure, gje lVtL / c ; and sim- 
ilarly for the kj field propagating in the other direction. 

The standard junction conditions at the front mirror 
imply that 

jj = VTf 3 + y/Rk 3 , g 3 = -VRfj + VTkj . (B8) 

We denote by X(t) the change of arm length produced 
by radiation pressure and the gravitational waves, and 
by X its Fourier transform. The oscillating X(t) pumps 
carrier light into the side bands. More specifically, in 
traveling from the front mirror z = to the perturbed 
position z = L + X(t) of the back mirror, then reflecting 
and propagating to the unperturbed location z = L, the 



carrier field acquires the formf^| 
^2 




\/2Dcos (u> [t - 2X(t)/c]) 



A_kTiuj A ( 2u u 

— I, — \ 7^D cosuj t + smuj t 

Ac V I \ c 



2tt J 
(B9) 



Comparing with the standard expression for th e fie ld at 
the location of the unperturbed end mirror [Eq. ( B2 ) with 
the amplitude changes indicated in the lower right of Fig. 
||, D -> (2/V2)(D/V2) and dj -> k j e-* lL l% we obtain 
the following expression for the field fed from the carrier 
D into the sideband amplitudes kj\ 



2 2uo n 
8ki=0, Sk 2 = — =D — -X . 

" Vt c 



(BIO) 



This acts as a source term in the standard junction con- 
dition for the back mirror: 



fcj-e 



-iQL/c 
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Tqj + 8k j 



(Bll) 



Note that qj is the noise-producing vacuum fluctuation 
that leaks into the cavity as a result of the optical losses. 



3. Cavity's internal field and radiation-pressure 
fluctuations 

By combining the front- mirr or and back-mirror junc- 
tion conditions (B8) and (Bll) we obtain for the side- 
band amplitude in the cavity 



VTfj + jRe mL ' c (Vf qj + 8k A 



(B12) 



Equations (g7j) and \fR = y/T^T = y/l - 0.033 ~ 1 
allow us to make the approximations \/Re inL / c ~ 1 in 
the numerator and [using (O)] 



(B13) 



1 - VRRe 2inL/c ~ (2L/c)(7„ - iQ) 



7 * = 7 (l + e/2) 



(B14) 



16 Here we have neglected the attenuation of the carrier field 
due to the arm-cavity losses. This neglect is in the same 
spirit as our ignoring attenuation in the input optics, in the 
beam splitter, and in mode matching into the arm cavities. 
Including these attenuations would simply change D in Eq. 



(B9) to D X (1 — ^power attenuation factor) 



Dx(l 



for the effect of arm-cavity losses. Equivalently, it would 
dictate replacing I by 7 (l— power attenuation factor) in IC, 
K,*, and all our formulas for the gravitational- wave noise. 
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in the denominator (accurate to better than 1 per cen t 
for all f2 of interest to us), thereby bringing Eq. (B12) 
into the form 



Jo 



(2i/c)(7,-tn) 



(B15) 



where we have used f = \tT . T he c avity's inter- 



nal electric field E- m t is expression (B2) with D — > 
(2/VT)(D/V2) [Eq. (p9h] and d 3 — ► jj [expression 
( |B15| )1; cf. Fig. |l6|. The power circulating in the cavity 
is this (Ef nt /4:ir)Ac, and consists of two parts, a steady 
classical piece 



w 14D \ 2 T 

Wch-c = --jrTiujo = —I 



Jo/2 
lL/c ' 



and a fluctuating piece 
6W ciTC -- 



y/I huj (fi + ^fe/2 qi) 



an 



(L/c)(7* 



2^ 



(B16) 

HC , 
(B17) 



where HC means Hermitian conjugate (adjoint) of the 
previous term. 



4. Mirror motion 



The circulating-power fluctuations (B17) produce a 
fluctuating radiation-pressure (back-action) force 



F BA = 2SW ciTC /c 



(B18) 



on each mirror. This force is equal and opposite on the 
cavity's two mirrors and, along with the gravitational 
waves, it produces the following acceleration of the mirror 
separation: 



d 2 X{t) 
dt 2 



1_ T d 2 h(t) 4SW cilc (t) 
dt 2 mc 



(B19) 



Here h(t) is the gravitational-wave field (projected onto 
the interferometer's arms), and rj ne is +1 for the north 
arm and —1 for the east arm (one arm is stretched while 
the other is squeezed). 

Below we will need an expression for the (Fourier trans- 
form of the) arm- length difference, x = X n — X e . It can 
be obtained by Fourier transforming the equation of mo- 
tion (B19 ), so lving for X (i.e., X n or X e ), inserting ex- 
pression"^ B17) for SW c j IC , and then taking the difference 
of the north and east arms. The result is 



x = Lh + xba 
[cf. Eq. dH)], where 

-4^/2I huj {ax + y/7j2 m) 



(B20) 



r7if2 2 L(7» — iCl) 



= -y/K*/2 LhsQh(a-L + \ftj2 m] 



J/3, 



(B21) 



Here we have introduced the quadrature amplitude for 
the difference of the arms' noise fields 

n n n e 

(B22) 



V2 



and have used Eq. @ for /f and /f, and Eqs. @, @, 
dl9| ) and ( |20| ) for the coupling constant /C*, the phase p*, 
the SQL power igQL and the standard quantum limit 

^SQL- 

Below we shall also need the following expression for 
the difference of the two arms' sideband fields produced 
by the mirror motions' coupling to the carrier: 



<5fc£ - 5k e 2 
~V2 



2y/2 



T V TiWo c 



(B23) 



This follows from Eqs. (B1C), (B3) and (ff 



5. Cavity output 

The field exiting from th e (n orth or east) cavity is ob- 
tained by combining Eqs. flBSj), (Bll) and (B12): 



\/Re 2iQL / c - VR , 
9j = r^r fj + 



1 - VRRe 2mL / c 



l TTq ] + VTSkj) e lUL/c 



1 - VRRe 2lQL / c 



(B24) 



Inserting Eq. ( B13| ) for the denominator and analogous 
expressions for the numerator, and discarding terms that 
are h igher order than linear in the losses, we bring Eq. 
(B24) into the form 



1 _ I A + ^g e i0 



9j 



where /?* is given by Eq. (J9q). 




(B25) 



6. Beam splitter output 



By combining Eqs. (|B4j), (JB25J) , and flB22| ), we obtain 
for the dark-port output of the beam splitter 

b 3 = ( 1 - h) a 3 e 2113 + VSme i0 



{c/2LfT ( Skf-5k] \ 
V ll + " 2 V V2 " 



(B26) 



Inserting <5fc™' e = [Eq. ( [B10| )1 and our expression (B22) 
for the difference of t he flfe 'S; and inserting Eqs. (B2C) 
for x and (|99|), (20), ( B14 ) for /C*, Tisql, 7*, we obtain 
for the output fields: 
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1 _ Is] aie W + y/£ nie if> , 



b 2 =(l- ^ a 2 e 2lf) + V£n 2 e^ 

"SQL 



(B27a) 



(B27b) 



By inserting expression (B21) for the back-action- 
induced mirror displacement xba, we obtain the input- 
output relations quoted in the text: Eqs. ( |97| ) and (101) 
with losses, and Eqs. (jig) in the lossless limit 



APPENDIX C: FILTER PARAMETERS 

In our discussion of FD homodyne detection [Sec. [v), 
we derived the following requirement for the conventional 
homodyne phase 9 and the filter parameters £j and Sj 
(with J — I and II): 



tan $ (ft) 



ft 2 ( 7 2 + ft 2 



A 4 



= tan 19 - 



an 



ai_ 



am 



aii_ 



(CI) 



[Eqs. ©, Q, and ©], where 

a J± = arctan(0 ± fl/Sj) (C2) 

[Eq. d§^)]. In this appendix, we shall show that this 
requirement is satisfied by the parameter choices asserted 
in the text: Eqs. (f|) and ©. 

We initially regard the par ameters 9,fj and Sj as un- 
known. By inserting Eq. ( |C2| ) into Eq. ( |Cl| ) and invoking 
some trigonometric identities, we obtain the requirement 

(Ro - Io cot 9) + (R 2 - h cot 6>)ft 2 + R A ft 4 



(i?o cot 9 + I a ) + (R 2 cot 9 + I 2 )9? + R 4 cot 9 ft 4 
7 2 ft 2 + ft 4 



A 4 



(C3) 



Here Rq + R 2 Sl 2 + i? 4 ft 4 is the real part and Iq + h^i 2 
is the imaginary part of (1 + i tan «i+)(l + i tanai_)(l + 
itanan+)(l +itanaii_). More specifically, 



Ro = i - g - & - mu + 

R2 = (l-g)/6l 1 + (l-&)/6 2 , 
R, = l/tffa) , 

h = 2Zu/6 2 + 2^/S 2 , . 



(C4a) 
(C4b) 
(C4c) 
(C4d) 
(C4e) 



To get rid of the ft 4 term in the denominator of Eq. 
(|C3|), we must set 



tt/2 , socot6l = 



(C5) 



(We cannot set i?4 = since that would require an infi- 
nite bandwidth for one or both of the filters.) To get rid 



of the ft 2 term in the denominator and the constant term 
in the numerator, and to make the ft 2 and ft 4 terms in 
the numerator have the correct coefficients, we must set 



J 2 = 0, 
Ro = 0, 
R 2 /(I R 4 ) = 7 4 /A 4 = 4/P 



R 2 j Ri 



7 



(C6a) 
(C6b) 
(C6c) 
(C6d) 



Here we have use d definition (89a) of the constant P. 

Equations (C6) are four equations for the four unkown 
filter parameters: the fractional frequency offsets £i, £n 
and the half bandwidths Si, Sn. In the next four para- 
graphs we shall explore the consequences of these four 
equations, arriving finally at the solution (89) for £j, (n, 
Si and Su gi ven in the text. 

Equation (C6a) implies that 



S?/S?i 



-W6 



(C7) 



Equation ( |C6b| ) implies that (1 - 66i) 2 = (£i +6i) 2 - 
It turns out that one of the frequency offsets is positive 
and the other is negative (cf. Fig. [l(]); we choose £i to be 
the positive one. It also turns out that £i + £n is positive 
(cf. Fig. [hJ). Consequently, we can take the square root 
of the above equation to obtain 



l - 66i = 6 + £n , 



(C8) 



which enables us to express the frequency offsets in terms 
of each other: 



6 = 



i + 6i 



1-6 

1 + 6 



(C9) 



Equation ( |C6c| ), when combined with Eqs. (07) and 
implies that 



(6 + 61) 2 



(CIO) 



We shall now combine this equation with Eqs. ( |C9| ) to 
obtain Eqs. ( |39| ) for the frequency offsets 6 and £n m 
terms of P — 4 j 4 / A 4 . Our first step is to define A± by 
Eqs. ( |89cj ) and flS9d| ), which are equivalent to 



.4, 



6i 



6 2 i - 1 



(Oil) 



Note that the relation (09) between £i and £n is equiva- 
lent to 



4A, A. 



1 . 



(C12) 



By using Eqs. (|C^), ( C 1 1 j ) and ( |C12| ), we can reexpress 
the right side of Eq. ( C10Q solely in terms of A + : 



8_ _ (AA\ - l) 2 
P ~ A+(4A 2 + + 1) 



(C13) 
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It is convenient to define Q by Eqs. (89b), which are 
equivalent to 



A, 



A. 



2Q/P 



(C14) 



Using Eqs. (|C12j) and (C14), we can rewrite (C13) in 
terms of Q instead of A + : 



2 
P 



2Q 
P 



P 

2Q ' 



which can be solved for Q as a function of P 

1 + yi + P 2 
q- 2 

This is the relation asserted in the text, Eq 



it completes our derivation of Eqs. (B9a)-( 39d| ) for the 
frequency offsets £i and £n in terms of P. 



Turn, finally, to the consequences of Eq. (C6d), which 
says 



By eliminating 5\\ with the aid of Eq. (| 



7 2 = ^i 2 6 



6 



(C17) 
we obtain 

(C18) 



Using Eqs. ( Cll ), ( C12 ), and ( 89b ), we can rewrite this 
as 



Si 



P 



(C19) 



which is the f ormu la for the half bandwidth Si given in 
the t ext , Eq. ( |S9cj ) . The correspond ing f ormul a fo r 5u, 
Eq. (§9$, follows directly from Eqs. jC19|) and @. 
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